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INTRODUCTION  AND  SUMMARY 


There  are  essentially  only  two  hypotheses  for  explaining  the 
narrow  V-wakes  of  ships  seen  by  SEASAT.  One  is  that  the  wakes  are  due 
to  straining  of  ambient  surface  waves  by  the  internal  wave  wake  of  the 
ship.  The  other  is  that  the  V-wakes  are  produced  by  the  Kelvin  wake 
of  the  ship.  (Here  the  word  Kelvin  wake  is  used  in  the  most  general 
sense,  and  includes  non-linear  processes  in  the  wake  as  well  as 
possible  interactions  with  the  turbulent  wake.  It  seems  clear  from 
the  data  that  such  non-linear  effects  are  important  in  the  near  axis 
Kelvin  wake). 

As  a  result  of  the  Georgia  Strait  experiment  in  the  summer  of 
1983,  and  the  Gulf  of  Alaska  experiment  in  the  spring  of  1984,  there 
is  now  a  considerable  amount  of  important  new  data  available.  The 
exercises  at  Georgia  Strait  consisted  of  two  separate  experiments. 

One,  conducted  ac  Knight  Inlet  where  the  buoyancy  frequency  is 
extremely  high  and  the  mixed  layer  is  shallow,  was  explicitly  designed 
to  produce  an  internal  wave  wake.  The  observed  (internal  wave)  wake 
was  consistent  with  calculations  based  on  standard  internal  wave 
theory  and  the  measured  density  profile.  The  other  experiment  was 
done  at  Dabob  Bay  under  conditions  that  are  closer  to  those  encoun¬ 
tered  in  the  open  ocean.  At  Dabob  Bay,  a  measureable  internal  wave 


wake  was  neither  expected  nor  observed.  However  the  SAR  did  see 
narrow  V  -  wakes.  Direct  observation  of  the  the  ocean  surface  clearly 
showed  that  these  narrow  V's  were  created  by  surface  waves  coming 
directly  from  the  ship. 

The  Gulf  of  Alaska  experiment  took  place  in  the  open  ocean  in 
a  region  with  an  extremely  deep  mixed  layer.  The  buoyancy  frequency 
was  essentially  zero  down  to  100m,  at  which  depth  there  was  a  peak  of 
about  9  cycles  per  hour.  No  visible  internal  wave  wake  should  be 
expected  under  those  conditions,  according  to  standard  theory.  Yet 
narrow  V-wakes  were  again  observed  by  the  SAR. 

The  experimental  results  appear  to  rule  out  the  internal  wave 
hypothesis  as  an  explanation  of  the  SEASAT  wakes.  They  also  seem  to 
indicate  that  the  actual  Kelvin  wake,  in  the  region  of  the  SAR  return, 
is  rather  different  from  what  the  conventional  Kelvin  wake  theory 
would  predict.  This  means  that  no  good  theory  exists  with  which  to 
explain  the  observed  results  quantitatively,  both  in  the  experiments 
and  in  the  SEASAT  images.  Nevertheless,  guided  by  the  experimental 
observations,  a  model  which  seems  at  least  qualitatively  reasonable 
can  be  constructed.  It  is  described  in  Section  III  of  this  report. 


The  Dabob  Bay  experiment  observed  narrow  V-wakes  during  a 
period  when  the  wind  speed  was  low  and  there  was  little  ambient  wave 
activity.  This  is  a  special  circumstance  for  two  reasons:  first  the 
pattern  in  the  Kelvin  wake  lasts  longer  and  secondly  the  SAR  directly 
sees  the  slopes  in  the  Kelvin  wake  rather  than  their  indirect  effect 
on  ambient  ripples.  Many  (and  perhaps  even  all)  bright  narrow  wakes 
observed  by  SEASAT  were  also  in  calm  seas.  Our  working  hypothesis  is 
that  Kelvin  wakes  appear  as  narrow  V's  only  when  the  wind  speed  is 
low.  This  hypothesis,  which  is  strongly  suggested  by  theory,  receives 
further  support  from  the  Gulf  of  Alaska  experiment.  There,  narrow 
V-wakes  seem  to  disappear  into  the  background  as  the  wind  speed 
increases . 

The  theoretical  model  described  later  in  this  report  is 
applicable  only  in  low  wind  speed  conditions.  Further  work  is 
required  to  understand  the  wind  dependence;  in  particular  the 
lifetimes  of  surface  waves  in  varying  wind  conditions  needs  to  be 
understood  better.  Some  progress  toward  this  goal  is  described  later 
(see  Appendix  C). 

The  generation  of  ship  wakes  (surface  or  internal)  and  their 
detection  by  radar  is  a  complicated  phenomenology  involving  various 
branches  of  physics,  hydrodynamics  and  oceanography.  Nevertheless, 
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with  one  notable  exception  (discussed  below),  we  believe  that  the 
theoretical  understanding  of  the  problem  is  adequate.  In  most  cases 
careful  calculations  should  be  accurate  to  within  a  factor  of  two. 

The  exception  is  the  inner,  small  angle  part  of  the  Kelvin 
wake.  In  this  region  the  wake  consists  of  short  wavelength  (typically 
tens  of  centimeters)  waves  travelling  almost  directly  outward.  At 
present  there  is  no  useful  theoretical  method  for  predicting  the 
amplitude  of  these  waves  or  even  their  gross  pattern.  Theoretical 
uncertainties  are  at  least  by  an  order  of  magnitude.  This  is  in 
contrast  to  the  situation  for  the  outer  part  of  the  Kelvin  wake  near 
the  cusp  at  19.5°  (half  angle).  Here  it  is  possible  to  make  rather 
accurate  theoretical  calculations.  Essentially  all  the  wave  drag 
comes  from  the  outer  part  of  the  wake. 

Any  theoretical  attack,  either  analytic  or  numeric,  on  the 
small  angle  Kelvin  wake  is  bound  to  be  a  difficult  undertaking. 
However,  the  problem  is  a  relative  „  unexplored  one  and  research  by 
qualified  investigators  should  be  encouraged. 

For  the  record,  the  remainder  of  the  theory  (which  is  in  good 
shape)  includes  (i)  the  generation  of  long  wave  length  (X  >  100m) 
internal  waves  by  moving  bodies  (short  wave  length  internal  waves  do 
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not  penetrate  the  mixed  layer  and  are  irrelevant),  (ii)  modulation  of 
surface  ripples  by  internal  waves  or  long  surface  waves  and  (iii)  the 
general  theory  of  radar  scattering  from  the  ocean  surface,  in  the 
regime  where  the  angle  of  incidence  of  the  radar  is  not  too  near 
grazing  and  the  modulation  of  surface  ripples  is  relatively  uaak. 

This  is  the  regime  of  interest  for  SEASAT  and  the  other  experiments 
discussed  here.  Under  these  conditions  no  measurements  incompatible 
with  the  theory  are  known  to  us. 

Part  I  of  this  report  briefly  reviews  the  experimental 
situation.  Part  II  discusses  the  evidence  regarding  the  Internal  Wave 
hypothesis.  Part  III  discusses  the  Kelvin  Wake  Hypothesis.  It 
describes,  first,  the  Dabob  Bay  results  and  a  theoretical  generalized 
Kelvin  wake  model  suggested  by  them;  and,  c.econd ,  some  of  the  SEASAT 
and  Gulf  of  Alaska  observations  in  the  light  of  the  Kelvin  wake 
hypothesis.  Part  IV  suggests  some  future  experiments. 

There  are  also  several  appendices  providing  further  details 


on  various  topics  alluded  to  in  the  text 


1.0  SUMMARY  OP  THE  EXPERIMENTAL  SITUATION 

The  narrow  V  wakes  which  were  first  observed  behind  surface 
ships  on  SEASAT  images  have  now  been  reproduced  under  much  more 
controlled  conditions  in  several  experiments.  There  appears  to  be  a 
remarkable  uniformity  among  the  results,  (except  for  Knight  Inlet, 
described  below)  which  can  be  summarized  by  the  following  list  of 
characteristics . 

1.  The  V-wakes  appear  as  narrow  bright  lines  behind  the 
ship.  Figures  1-la,  1-lb,  and  1-lc  give  examples  from 
SEASAT  (Rev.  407),  Dabob  Bay  (DREP  8/3),  and  the  Gulf  of 
Alaska  (Ship  pass  006,  03/13/84).  The  similarity  is 
striking. 

2*  In  all  cases  where  the  ship  speed  is  known,  the  V-wake 
half  angles  6  (after  various  corrections  are  made)  are 
consistent  with  the  formula  tan  6  »  +  C  /U  cos  \J»  , 

O 

where  C  *8  t'ie  8ur^ace  wave  g^oup  spaced  appropriate  to 
g 

the  SAR  Bragg  wavelength,  is  the  angle  between  the 
ship's  course  and  the  SAR  flight  direction,  and  U  is  the 
ship  speed. 


Figure  1-lb.  V-wake  in  Dabob  Bay. 
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The  pattern  of  radar  cross  section  change  across  the 
wakes  is  also  very  similar,  as  shown  in  Figures  l-2a  and 
2b  (again  for  407,  DREP  8/3);  the  results  in  the  Gulf  of 
Alaska  look  similar.  There  are  plots  of 
Ao/o  (actually  /Ao/o  for  407)  across  the  wake  about  1 
km  behind  the  ship .  The  peaks  in  A o/ a  are  large 
(Ao/o  is  about  10  in  Dabob  Bay  and  4  in  407)  and 
generally  quite  narrow  (at  most  a  few  pixels). 

The  wake  lengths  vary,  up  to  perhaps  20  km,  and  most  if 
not  all  seem  to  occur  at  low  wind  speeds.  In  both  Dabob 
Bay  and  the  Gulf  of  Alaska  wind  speeds  were  known,  and 
the  V's  disappeared  for  wind  speeds  above  2-3  m/sec.  In 
SEASAT  wind  speeds  were  not  known,  but,  judging  by  the 
dark  backgrounds  where  tracks  are  seen,  were  apparently 
low  when  wakes  were  visible. 

The  wakes  tend  to  disappear  as  the  angle  <(>  between  the 
SAR  track  and  the  ship  track  Increases.  Also,  the  V- 
wake  angle  9  seems,  from  Dabob  Bay  data,  to  decrease  as 
iji  increases. 


The  wakes  (shown  by  Doppler  offset  in  cases  where  ij;  is 
large)  are  approximately  stationary  in  the  water  rest 
frame,  not  in  the  ship  rest  frame:  they  do  not  travel 
with  the  ship. 

At  least  some  wakes  show  multiple  crests  (SGASAT  Rev. 

407  and  perhaps  some  other  SEASAT  wakes;  also, 
apparently,  Gulf  of  Alaska  wakes.  (There  is  an 
indication  of  other  peaks  in  the  plot  of  A o/o  vs.  0  for 
Dabob  Bay) . 

One  case  (Knight  Inlet,  DREP  10/8)  looks  rather 
different  from  the  others.  There  are  many  crests,  and 
the  cross  section  modulation  Ao/o  is  very  small.  Both 
the  V  angle  and  the  width  of  the  pattern  are  much  larger 
in  the  other  cases.  DREP  10/8,  we  remind  the  reader, 
was  a  slow  ship  (5  knots)  moving  in  a  very  highly 
stratified  region  (peak  buoyancy  frequency  of  69  cph  at 
5  m  depth),  for  the  deliberate  purpose  of  exciting  an 
internal  wave  wake. 

The  SAR  look  angle  i  from  vertical  appears  to  be  a 
crucial  parameter.  SEASAT  was  at  23°  off  vertical,  and 
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most  of  the  runs  at  Dabob  Bay  and  in  the  Gulf  of  Alaska 
attempted  to  stay  close  to  this.  In  one  case,  however, 
in  the  Gulf  of  Alaska  i  was  changed  from  about  30° , 
where  the  V  wake  was  plainly  visible,  to  about  50°, 
where  the  same  wake  entirely  disappeared.  It  is  worth 
noting  that  SIR  A,  a  SAR  flown  on  the  space  shuttle,  saw 
surface  ships  but  no  narrow  V  wakes.  It's  look  angle 
was  also  above  50°. 

10.  If  the  radar  cross  section  is  due  to  Bragg  scattering, 
(at  the  SEASAT  angle  of  23°  and  at  the  lower  angles  in 
Dabob  Bay  and  the  Gulf  of  Alaska)  the  scattering  is  from 
surface  waves  of  around  30  cm  wavelength.  In  Dabob  Bay 
the  slope  spectra  across  the  wake  were  measured,  and 
optical  photographs  of  the  wake  were  also  taken. 

Neither  of  these  measurements  seems  to  show  much 
enhancement  over  background  at  30  cm  wavelength  (though 
both  show  considerable  enhancement  at  larger 
wavelengths) . 
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2.0  INTERNAL  WAVE  HYPOTHESIS 


The  hypothesis  that  ship  generated  internal  waves  are 
responsible  for  the  narrow  V-wakes  observed  behind  surface  ships  by 
SEASAT  appears  to  be  incompatible  with  existing  experimental 
results.  This  conclusion  is  based  on  the  following  observations. 

A.  Conventional  theory  predicts  that  the  radar  signature 
from  straining  of  ambient  surface  waves  by  ship  generated  internal 
waves  is  too  small  by  two  or  three  orders  of  magnitude  to  be  visible 
on  SEASAT. 


There  are  now  several  detailed  numerical  calculations  of 
radar  cross  section  modulations  for  several  different  pycnoclines  and 
surface  ships,  appropriate  to  SEASAT  Rev.  407,  to  Dabob  Bay,  to 
Knight  Inlet,  and  to  the  Gulf  of  Alaska.  All  of  these  are  pretty 
well  consistent  among  themselves,  and  also  agree  reasonably  well  with 
earlier  simpler  theoretical  calculations.  Predictions  of  these 
calculations  for  various  relevant  quantities  in  the  three  situations 
are  listed  in  Table  2-1. 

B.  Conventional  theory  is  validated  by  the  following 
sequence  of  experimental  results. 
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TABLE  2-1 


EXPT.  N  (PEAK)  DEPTH  OF  PEAK  WATER  DEPTH 


DREP  8/3 

.06 

N 

O 

— 

9/2 

.06 

.03 

— 

- 

10/8 

— 

— 

1.4 

.4 

Gulf  of  Alaska 
(Ship  Pass  006, 

3/13/84) 

— 

.005 

- 

NOTE ;  The  three 

407  numbers  are 

three  different 

guesses  as 

to  the 

(unmeasured)  profile. 


a.  As  predicted  by  the  theory,  an  internal  wave  wake  was 
observed  behind  both  QUAPAW  and  ENDEAVOUR  in  Knight  Inlet  (DREP 
10/8).  This  wake  is  quantitatively  reasonably  consistent  with  the 
theoretical  values  for  radar  return  and  surface  currents,  as 
indicated  in  Table  2-1.  Further,  the  observed  geometrical  pattern  of 
the  wake  closely  resembles  the  predicted  pattern. 

b.  As  predicted  by  the  theory,  no  internal  wave  wake  was 
observed  in  Dabob  Bay.  That  the  narrow  wake  seen  in  Dabob  Bay  is  not 
an  internal  wave  wake  is  indicated  by  the  following  observations. 

1.  The  surface  wave  slopes,  as  measured  directly  from 
ENDEAVOUR,  show  an  enhancement  around  the  Kelvin  wake 
wavelength  as  a  function  of  wake  angle  from  19  V2 0  down  to 
about  2V20>  at  which  point  everything  disappears, 
presumably  due  to  wave-wave  interactions  among  the  very 
short  waves,  or  to  interactions  with  the  turbulent  wake.  At 
the  angle  at  which  the  V-wake  is  seen  (2.6°),  there  is  no 
'  evident  enhancement  over  a  wide  band  of  surface  wave 
wavelengths,  such  as  one  would  expect  to  be  produced  by 
internal  wave  straining.  Only  the  enhancement  near  the 
Kelvin  wake  wavelength  is  seen. 


V-wake  in  Dabob  Bay  show  very  large  (factor  of  10) 
enhancements  at  +2.6°  from  the  wake  axis.  These 


is  far  larger  than  the  observed  upper  limit  of  22  m  for  the 
signal  width.  Couversely,  for  an  internal  wave  of  100  m 
wavelength,  the  group  speed  is  only  20  cm/sec,  and  the  phase 
speed  is  only  34  cm/ sec.  These  are  much  too  slow  to  produce 
a  signal  at  2.6°» 

For  an  internal  wave  to  produce  an  effect  over  a  more 
localized  region  than  a  quarter  wave  length  requires  that 
one  invoke  a  more  complicated  surface  wave-internal  wave 
interaction  mechanism,  for  which  there  seems  to  be  no  other 
evidence. 

c.  Much  experimental  information  on  natural  internal  waves 
exists.  None  of  it  appears  to  be  incompatible  with  the 
standard  theory,  as  described  below. 

One  of  the  surprises  of  the  SEASAT  SAR  experiment  was 
that  internal  wave  images  were  observed  frequently  along 
coasts,  over  sea  mounts  and  even  in  the  open  ocean.  This 
surprise  stems  from  the  fact  that  internal  waves  are 
essentially  motions  internal  to  the  sea,  having  maximum 
particle  displacements  where  the  maximum  vertical  density 
gradients  are  found,  10's  to  100's  of  meters  below  the 
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surface.  Nevertheless,  surface  manifestations  of  Internal 
waves  are  well  observed  by  spacecraft  and  aircraft 
photography  as  well  as  SAR  Images.  The  large  scale  and 
unique  shape  of  internal  wave  packets  and  the  equally  large 
scale  SEASAT  SAR  images  (100  km  wide  swath  some  100's  to 
1000' s  of  km  long)  account  in  some  measure  for  the  frequent 
observations. 

The  permanence  of  internal  wave  packets  in  SAR  images 
has  prompted  a  number  .of  hypotheses  to  explain  how  Internal 
waves  couple  to  the  surface  roughness  to  which  the  radar 
waves  are  sensitive.  Below  we  list  several  of  these 
hypotheses: 

1.  Enhanced  surface  roughness  in  current  convergence 
zones. 

2.  Surface  current  modulation  of  surfactants. 

3.  Surface  current  modulation  of  SAR  image  intensity  via 
doppler  effects. 
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4.  Surface  current  modulation  of  wind  generation  of  small 
surface  waves  (  ~30  cm  wavelength). 

Quantitative  theory  for  these  mechanisms  is  available  for 
items  1  and  3  and  we  shall  emphasize  them  in  our  discussion 
below. 


Several  SEASAT  and  aircraft  SAR  images  of  natural 
internal  waves  have  been  analyzed.  In  some  cases 
corresponding  surface  and  sub-surface  measurements  have  also 
been  made.  Below  we  list  several  "experiments"  in  which  SAR 
images  of  internal  waves  have  been  analyzed,  noting 
location,  reference  and  extent  of  surface  measurements. 

1.  Georgia  Straits  (1978  &  1983):  Aircraft  and  satellite 
SAR  observations  of  internal  waves  plus  extensive 
surface  measurements  were  done  in  1978  with  further 
similar  experiments  involving  aircraft  SAR  in  1983. 
(Hughes  &  Grant,  1978;  Hughes,  1978  and  Hughes  &  Gover, 
1983). 

2.  Massachusetts  Bay  (1978):  Four  SEASAT  SAR  passes  were 
analyzed  in  conjunction  with  buoy  data  from  near 
Stellwagen  Bank  (Trask  and  Briscol,  1983). 
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3.  Gulf  of  California:  Nine  SEASAT  SAR  images  were  used 
to  study  internal  waves  generated  by  tides.  Simple, 
non-linear  theory  allows  estimates  of  internal  wave 
amplitude  and  tidal  losses.  No  surface  observations 
were  made  (Fu  &  Holt,  1984). 

4.  JASIN  (1978):  Seventeen  SEASAT  SAR  passes  were  made 
over  the  Joint  Air-Sea  Interaction  (JASIN)  experiment 
area  west  of  Scotland  in  summary  1978.  Internal  waves 
were  observed  on  numerous  occassions  and  classified  in 
terms  of  location,  local  wind  speed,  water  depth, 

etc.  Considerable  surface  and  subsurface  data  was  also 
collected.  However,  the  SAR  and  in  situ  observations 
have  not  yet  been  analyzed  together  (Kasischke,  et  al, 
1983;  Levine  et  al.,  1983  and  Bagg,  1983). 

We  point  out  here  that  research  concerning  the 
observation  of  shallow  water  bathymetric  features  by  SAR  are 
relevant  to  the  internal  wave  problem.  Much  of  the  same 
surface  current-surface  roughness  coupling  physics  probably 
occurs  in  both  cases.  The  source  of  surface  currents  is,  of 
course,  different  in  the  two  cases.  Some  recent  papers  are 
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published  in  the  Proceedings  of  the  Conference  on  the  Future 
of  Remote  Sensing  of  the  Troposphere  and  Oceans  held  in 
Israel  in  May  1984. 

The  primary  question  we  ask  of  these  data  sets  is,  are 
they  consistent  with  current  understanding  of  the  mechanism 
by  which  SAR  can  image  surface  manifestations  of  internal 
waves.  By  current  understanding  we  refer  to  item  1  in  the 
list  of  possible  mechanisms  above,  namely  that  the  internal 
wave  induced  surface  currents  interact  with  short  (  '•*30  cm  ) 
surface  waves  via  weak  hydrodynamic  interaction.  This  leads 
to  modulation  of  radar  cross  section  (a)  by  typical 
internal  waves  of  (6o/o)  on  the  order  of  unity.  Although 
the  data  sets  so  far  analyzed  are  not  comprehensive  (about  a 
dozen  cases  quantitative  comparisons  have  been  made),  they 
are  generally  in  reasonable  agreement  with  the  above  theory. 
There  were  observations  of  natural  internal  wave  features 
during  the  1983  Georgia  Strait  experiment  which  showed 
(6o/o)  "•  10  .  However,  these  cases  (one  in  particular 
called  the  "hum  dinger")  were  exceptionally  strong,  internal 
bores,  expected  to  be  strongly  non-linear. 


d.  Further  support  for  the  standard  theory  comes  from  an 
analysis  of  a  Kelvin  wake  in  a  SEASAT  image  processed  by  RAE 
Farnborough,  shown  in  Figure  2-1.  The  Farnborough  picture 
shows  a  range-travelling  ship  plus  an  associated  39°  (full 
angle)  wake.  The  wake  shows  a  pattern  of  transverse  waves 
filling  the  interior  of  the  39°  wedge,  rather  than  the  two 
very  narrow  bright  rays  visible  on  most  SEASAT  pictures. 
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Figure  2-1.  The  Famborough  Picture. 
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The  observed  pattern  seems  to  fit  the  structure  of  the 
standard  Kelvin  wake,  suggesting  that  the  water  motions 
responsible  for  the  radar  cross-section  modulation  are 
nothing  more  or  less  complicated  than  the  Kelvin  wake.  (It 
is  interesting  to  note  that  the  observation  of  a  range 
travelling  ship  in  the  Gulf  of  Alaska  (ship  pass  010. 
3/13/84)  also  shows  the  usual  Kelvin  wake  but  not  the  bright 
narrow  V  r7ake.)  This  gives  us  an  opportunity  quantitatively 
to  test  the  hypothesis  that  the  modulation  is  due  to 
straining  of  an  isotropic  background  of  surface  waves  (or 
due  to  tilting  such  a  background  by  a  flow  pattern  produced 
by  the  ship). 

Let  us  first  establish  some  kinematic  and  geometrical 
facts  about  this  picture.  Horizontal  scale  is  38  m/mm  on 
the  picture.  The  transverse  wave  system  lias  about  4 
crests/lO  mm  or  an  apparent  wavelength  of  95  m.  The 
wavelength  of  the  transverse  system  behind  a  ship  travelling 
at  speed  U  is 


The  implied  ship  speed  is  U  =  12  m/s  »  24  knots  (1)  A  check 
on  this  comes  from  the  displacement  of  the  radar  image  of 
the  ship  from  the  head  of  the  V.  The  vertical  scale  of  the 
picture  is  37  m/mm.  The  ship  offset  is  15  mm  or 
Ax  *  550  m  .  The  in-range  component  of  the  ship's  velocity 
(from  the  picture,  in  fact,  the  ship  appears  to  be 
travelling  precisely  in  range)  is 


U  -  Ax 


SAT 


R  sin  0 


inc 


where  Vg^,,  R  are  the  satellite's  velocity  and  range.  This 
gives  U  *  12.7  m/s,  which  is  consistent  with  the  wavelength 
of  the  Kelvin  wake. 


Next,  we  want  to  estimate  the  amplitude  of  the  Kelvin 
wave  train  in  order  to  determine  the  strain  and  slope  field 
which  the  radar  is  seeing.  The  far  field  transverse  wave 
field  can  be  fairly  reliably  calculated  from  slender  ship 
theory  at  least  far  from  the  wake  axis.  The  answer  is,  for 
a  ship  whose  beam  to  length  ratio  is  1/10, 

„2 

displacement  5  ■  .6  —  /  —  cos  (^v) 


slope  m  =  .6  /  —  cos  (—■) 

x  x  uz 


W  T 

strain  rate  ^  —  cos  (^r) 

3x  U  x  y2 


where  x  is  the  down  track  distance  and  L  is  the  ship  length. 


The  radar  image  suggests  that  the  ship  length  is  about 
240  m.  For  x  »  6  km,  the  maximum  extent  of  the  Farnborough 
picture,  £  ~  .25  m  and  ~  1.  Since  slender  ship  theory 
probably  over  estimates  the  strength  of  the  wake,  these 
numbers  should  be  decreased  by  a  fudge  factor  of  a  few. 
Ideally,  the  wake  could  be  computed  using  Carl  Scragg's 
computer  code . 


Let  us  now  consider  the  strain  issue:  The  period  of  the 
transverse  wave  system  is 


2""2/8  V 

T  -  -.T—a-  °  2i^  -  7.7  sec 
U  g 


The  maximum  strain  in  half  a  period  is 


e  -  d)  <^)  -  (*r)  (-6  |  ) 


(.6  it)  ✓  - 
x 


-  .38  @  x  ■  6000  m 

The  strain  modifies  the  ambient  surface  wave  spectrum  and  we 
have  to  decide  whether  we  are  in  the  strain  or  strain  rate 
regime  (i.e.,  is  the  lifetime  of  the  relevant  surface  waves 
long  or  short  compared  to  the  period  of  the  Kelvin 
transverse  waves).  The  Bragg  wavelength  is  30  cm  and  the 
period  of  such  waves  is 

T  ■  /  —  X  ■  .44  sec 
g 

The  Kelvin  wave  period  is  17  times  as  long.  We  will  assume 
that  this  is  small  enough  that  decay  of  the  surface  wave  is 
negligible  and  we  are  in  the  strain  regime. 

The  empirical  formula  for  radar  backscatter  modulation 
by  a  strain  field  e(x)  ,  in  the  strain  dominated  regime,  is 

— ■  c  e  (x) 
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where  c  is  nearly  one.  The  strain  field  of  the  transverse 
Kelvin  wake,  according  to  the  preceding  paragraph,  is 

e(x)  -  .38  /  cos  (^) 

x  ^ 

It  is  perfectly  believable  that  the  pattern  in  the 
Farnborough  picture  is  due  to  a  40%  intensity  modulation. 

It  would  be  good  to  have  the  digitized  data  so  as  to  measure 
the  intensity  modulation. 

The  first  order  conclusion  is  that  the  picture  of  SAR 
returns  as  arising  from  strain  modulation  of  wind  generated 
ambient  waves  fits  quantitatively  what  seems  to  be  going  on 
in  the  Farnborough  picture.  The  key  point  is  that  in  this 
picture  it  is  clear  that  the  Kelvin  wake  is  directly 
responsible  for  the  radar  return  and  that  there  is  a  well 
developed  ambient  wave  spectrum  for  the  Kelvin  wake  to 
modulate.  It  would  obviously  be  of  great  interest  to 
combine  Scragg's  wake  code  with  a  radar  cross-section  model 
to  obtain  a  quantitative  prediction  of  the  radar  return  to 
compare  with  the  Farnborough  picture. 


2-16 


e.  To  explain  the  3EASAT  images  by  internal  wave  wakes 
requires  that  surface  ships  radiate  internal  waves  much  more 
efficiently  than  u-s  been  believed.  One  way  to  bound  such  a 
hypothetical  mechanism  is  to  compute  the  drag  that  would  be 
associated  with  the  internal  wake.  This  has  been  done  for 
selected  cases  with  the  results  shown  in  Table  2-2.  The 
power  required  to  overcome  IW  drag  (last  column)  is 
reasonable  in  all  cases  where  the  standard  theory  has  been 
used  (all  rows  except  the  bottom  two).  In  the  bottom  two 
rows,  the  amplitude  of  the  internal  waves  were  (arbitrarily) 
multiplied  by  a  factor  of  7  and  10  to  bring  the  radar  cross 
section  up  to  a  level  consistent  with  the  SEASAT 
observations,  giving  8,000  and  16,000  hp  of  IW  drag.  Since 
large  commercial  ships  typically  have  power  plants  with 
30,000  hp,  we  believe  that  IW  drags  of  this  magnitude  are 
ruled  out  by  experience.  The  possibility  of  bounding 
anomalous  IW  wakes  by  putting  limits  on  drag  should  be 


studied  in  more  detail 


TABLE  2-2 


B-V  PROFILE 

SHIP 

SPEED 

ENERGY* 

IN  INTERNAL 
WAVES 

POWER* 

TO  OVERCOME 
IW  DRAG 

DREP  9/2 

QUAPAW 

DABOB  BAY 

Measured 

16  kts 

4.8  j/m 

.05  hp 

DREP  10/8 

QUAPAW 

KNIGHT'S  INLET 

Measured 

5  kts 

1,924  j/m 

6.5  hp 

(#1) 

Historical 

16  kts 

825  j/m 

8.9  hp 

(#2) 

SEASAT 

Optimistic 
(For  IW's) 

16  kts 

15,000  j/m 

160  hp 

REV.  407 

CAPE  CANAVERAL 

Optimistic 
with  source 
increased  to 
get  6o  ■  .7 

0 

16  kts 

730,000  j/m 

7,800  hp 

Increased  to 

get  6o  -  1  16  kts 

1,500,000  j/m 

16,000 

hp 

0 

Table  Caption;  Internal  wake  energies  computed  by  S.  Bovchardt  (Dynamics 
Techology).  Even  with  the  "optimistic"  profile  the  internal  wave  amplitudes 
for  SEASAT  Rev.  AO )  are  too  small  to  explain  the  observations.  To  obtain  the 
last  two  rows  the  amplitudes  are  (arbitrarily)  multiplied  by  7  and  10  to  get 
6 o/o  of  .7  and  1,  respectively. 


*  Keeping  10  modes  and  wavelengths  greater  than  100m,  the  true  IV 
drag  is  greater. 
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3.0  KELVIN  WAKE  HYPOTHESIS 


This  is  the  second  serious  contender  for  an  explanation  of 
the  narrow  V's.  There  are  two  obvious  things  to  be  said: 

A  Kelvin  wake  with  its  characteristic  pattern  of  divergent 
and  transverse  waves  does  not  look  at  all  like  a  narrow  V.  But  it  is 
readily  shown  that  in  a  classical  Kelvin  wake  the  wavenumber  vector 
k  is  a  constant  for  constant  wake  angle  6  (see  Figure  3-1).  Thus  if 
the  SAR  is  k  -selective  (as  it  is  in  Bragg  scattering)  this  would 
indeed  lead  to  V-type  signatures. 

We  do  not  understand  the  Kelvin  wake  theory  for  the  interior 
wake  (small  0  )  even  one  hundred  years  after  publication  of  Lord 
Kelvin's  classical  paper.  In  the  past  the  emphasis  was  on  the  outer 
19.5°  cusps  which  are  important  to  the  problem  of  wave  drag. 

Nevertheless  we  will  try  to  interpret  the  Dubob  results  in 
terms  of  classical  Kelvin  wake  theory.  (The  attempt  turns  out  not  to 
be  successful!  but  it  gives  important  clues.)  The  present  assumption 
is  that  we  see  a  classical  Bragg  reflection  from  the  divergent  waves 
of  a  classical  Kelvin  wake.  (We  believe  this  suggestion  first  came 
from  NOSC.)  But  we  must  say  at  once  that  we  believe  that  the  SEASAT 
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two  0-values. 


phenomenology  is  quite  different  for  weak  and  strong  winds  (indeed  it 
may  not  exist  for  strong  winds),  and  for  near-azimuth  and  non-azimuth 
travelling  ships;  Dabob  Bay  is  representative  only  of  the  case  of 
calm  and  more  or  less  azimuthal  ships . 

1*  The  angles  are  defined  in  Figure  3-2.  SAR  looks  down  at 
an  angle  i  from  the  vertical.  A  ship  in  the  field  of  view  travels  at 
an  angle  to  the  right  of  the  satellite  course.  The  SAR 
coordinates  £,ri  and  the  ship's  coordinates  x,y  are  shown;  both  are 
useful.  Note  that  the  ship  is  travelling  in  the  negative  x 
direction.  The  look  angle  (normal  to  the  satellite  course)  is 
iji  =  0,  180°  for  an  azimuth  travelling  ship,  and  ip  =  Hr  90 0  for  a 
range  travelling  ship. 

2.  The  classical  Kelvin  solution  consists  of  a  series  of 
waves  n  =  1,  2,  .  .  .  which  are  stationary  in  a  reference  frame 
fixed  to  the  ship.  In  this  frame  the  water  velocity  is  every  where  U 
in  the  positive  x-direction.  A  wave  crest  y(x)  forms  an  angle 
tan<j>  =>  iy/dx  with  the  x-axis.  The  flow  component  normal  to  the 
crest  is  U  siu<j>  ,  and  the  condition  of  stationarity  is  that  this 
should  be  equal  and  opposite  to  the  wave  phase  velocity  C  (see 
insert,  Figure  3-3): 
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Note:  The  ship  moves  at  an  angle  y  to  tfiu  right  of  SEASAT.  Thus  ip  =  80*  corresponds  to  a  positive  rate-rango,  henoo  a 
negative  Doppler,  and  In  SEASAT  coordinates  {,  y  this  leads  to  a  relative  ship  displacement  In  tho  positive  ( 
direction, 


Figure  3-2.  Definition  plot. 


Note:  The  ship  Is  located  at  0,  moving  at  U  =  5  m/s;  It  occupied  positions  A  and  A’  at  times  t  =  280s  and  240s  ago, 
respectively.  At  these  times  Impulsively  generated  waves  spread  outwardly.  For  a  Bragg  resonance  at  i  =  2m  (for 
Illustration  In  this  figure)  the  group  velocity  Is  Cg  =  0.88  m/s,  and  the  distance  traveled  In  times  t  =  280s  is  Cgt  = 
247  m,  as  shown.  (Longer  waves  are  further  out,  shorter  waves  are  closer  In.)  The  radar  look  angle  Is  if,  and  so  the 
location  for  Bragg  scatter  Is  In  a  small  spot  centerf.u  on  B.  Similarly  for  point  A'  the  scattering  Is  from  B’.  These 
points  lie  along  a  straight  line  OB  forming  the  starboard  wing  of  the  V,  at  an  angle  8  relative  to  the  ship's  track.  A 
similar  construction  applies  to  the  port  wing.  It  will  be  seen  that  port  8  Is  smaller  (for  if  as  shown)  and  that  both  8’s 
decrease  as  if  approaches  90*.  The  figure  also  shows  the  stationary  Kelvin  crest  pattern  for  wave  crests  n  =  10,  20, 
30,  40, 106,  and  123.  In  the  special  case  that  if  =  ♦  (the  slope  dy/dx  of  the  stationary  waves),  the  Bragg  resonance 
occurs  at  point  G  and  the  waves  are  stationary  with  respect  to  the  ship's  coordinates.  At  that  point  the  cyllndrlcally 
spreading  waves  are  tangent  to  the  stationary  Kelvin  waves. 


Figure  3*3.  Geometric  construction  for  SAR  Vs. 


Loci  of  equal  <J>  lie  along  a  straight  line  which  forms  an  angle 
8  with  the  wake  axis  (marked  by  the  turbulent  wake),  with 


tan8  -  - tan^_  (2 

2  +  tan  $ 

For  small  angles  8  *  V2  ♦  • 

Each  crest  n  has  a  divergent  wave  branch  and  a  transverse 
wave  branch;  these  join  at  the  line  of  cusps  where 

8  »  arc  sin  (1/3)  ■  19.47°  (3' 

We  are  concerned  chiefly  with  the  divergent  waves.  The  wave  length 
X  is  relatively  long  near  the  cusp,  and  diminishes  towards  the 
Interior  (small  8  and  small  $)  as 

k  u  __ — SL_ —  t  \  m  2-n/]n  ,  (4] 

U  sin  $ 

The  phase  velocity  C  ■  /  g/k  ■  U  sin$  in  agreement  with  (1),  so 
that  the  relation  (4)  can  be  regarded  as  a  purely  geometric 
consequence  of  stationarity.  The  wave  number  components  are 


k  *  k  cos<j>  ,  k  »  k  sin$  .  w7 

x  y 

For  the  transverse  wave  at  the  axJ  s  we  set 
4  =  90^,  hence  k  =  g/U^,  snd  C  *  -U. 

3.  The  Bragg  condition  states  that  backscatter  is 
pronounced  when  the  ocean  wave  crests  are  lined  up  normal  to  the  look 
angle, 


4  =»  U;  (b) 

9BRAGG  * 

with  a  wave  number  (for  first  order  Bragg) • 

^•21^  sin  i  .  (?) 

Thus  Bragg  resonance  from  the  Kelvin  wake  requires  a  very  special 
combination  of  four  independent  experimental  parameters  (see  Appendix 
B  as  well):  the  radar  frequeucy  fR,  the  radar  inclination  i,  the 
(relative)  ship's  course  4  and  speed  U.  From  (4),  (6),  and  (7), 
this  is 


1 


I 


2k^  sin  i 


2  2 
U  sin  4 


where  kR  is  the  radar  wave  number.  In  general  this  condition  cannot 
be  satisfied;  in  particular  only  one  side  o£  the  classical  Kelvin 
wake  can  be  lined  up  properly  so  that  only  one  side  of  a  V-wake 
should  be  visible.  We  conclude  that  the  classical  Kelvin  wake  cannot 
explain  the  narrow  V-wakes. 

4.  Suppose  instead  that  the  wave  generation  is  variable. 

We  can  model  it  by  a  series  of  charges  dropped  from  the  stern  of  the 
moving  vessel  (not  necessarily  at  regular  intervals)  (see  also 
Appendix  B).  In  the  limit  of  zero  interval  between  equal  charges  we 
recover  the  situation  for  a  stationary  phase  solution.  For 
definiteness  Figure  6  has  been  drawn  for  U  “  5  m/sec,  and  we  now 
consider  the  situation  1.4  km  and  1.2  km  behind  the  ship,  two  points 
that  were  occupied  by  the  ship  280s  and  240s  ago,  respectively.  From 
each  of  those  two  points  (as  from  any  other  point)  we  have  radial 
waves  spreading  outward,  with  group  velocity  Cg{  for  deep  water 
surface  waves 

C  -  1/2  C  -  V2  /  g'X/'i*  (8) 

O 

The  long  waves  are  fast  and  out  in  front,  the  short  waves  are  slow 
and  nearer  their  point  of  generation.  The  circular  arcs  are  drawn 


for  an  assumed  Bragg  wavelength  \  -  2.00  m  which  leads  to  C  -  0.353 
U.  The  arc  to  the  left  has  a  shorter  radius,  having  been  generated 
more  recently.  The  conditions  for  Bragg  resonance  occur  in  the 
vicinity  of  the  arcs  at  the  points  B,  B'  where  the  arcs  are  at  right 
angle  to  the  look  direction  i}i  .  For  the  triangle  OAB  we  have  OA  » 
Ut,  AB  -  Cgt,  and  hence 

C  cos  (i|>-0)  «  U  sin8  , 

g 

or 

C 

-£■  cos  \|) 

tan0  -  +  - ^ -  (9) 

1  +  -f  sin  * 

The  two  lines  correspond  to  the  two  wings  of  the  V;  a  positive 
9  corresponds  to  the  starboard  wing,  a  negative  6  to  the  port  wing. 

The  radar  wave  length  and  inclination  determines 
cg  "  I  8*RADAR^8*n  ^  u  is  shiP  speed,  and  \J>  is  the  angle 
between  satellite  and  ship  course.  Equation  (9)  then  predicts  the  V- 
angles  (Figure  3-3).  We  note  a  few  special  cases,  designated  by 
C,D,E  along  the  circular  arc: 


E:  Azimuth  travelling  ship,  =  0,  tan6  =  C  /U 

§ 

D:  Maximum  value  of  0,  ij>  =  8  ,  sin0  =  C  /U 

8 

C:  Stationary  phase  point,  circular  arc  tangent  to 

stationary  phase  crest.  At  this  point  the  stationary  phase 

condition  (1)  holds,  sin<f>  =  sinij)  =  C/U  =  2C  /U  ,  and  when 

8 

this  is  substituted  in  (9)  we  recover  the  stationary  phase 
result  (2); 

4)  =  <p  ,  tan0  =  — 

2+tan  <j> 

Presumably  the.  strongest  return  comes  from  the  vicinity  of  C,  because 
it  is  here  that  sources  very  close  to  A  interfere  constructively. 

Equations  (9)  predict  some  conditions  which  can  be  checked 
with  SAR  images; 

1.  The  wing  angles  are  constant  and  nearly  equal  at 

0  •>  +  C  /U  for  small  to  moderate  values  of  U)  . 

-  8 

2.  As  \{i  +  45°  there  is  a  slight  (but  measurable)  wing 


asymmetry 


3.  For  increasing  *  the  V  becomes  increasingly  narrower. 


4.  In  SAR  images  (which  uses  Doppler  for  determining  image 
coordinates)  there  is  a  combination  of  effects  which 
are  discussed  in  the  next  section. 

5.  For  an  increasing  range-rate  the  Doppler  is  negative  and 
the  object  is  displaced  in  the  direction  of  objects  from  which  SAR  is 
receding:  that  is  towards  positive  5  •  Thus  the  £  offset  is 
proportional  to 

U  sin*  for  the  ship, 

+  C  cos*  for  the  S/P  wings  . 

Accordingly  wings  and  ships  are  offset  equally  if 

tan*  -  +  C/U  «  +  2  C  /D  for  S/P  . 

-  -  g 


This  is  shown  in  Figure  3-4 


Note:  The  Bragg  curves  are  drawn  for  a  ship  speed  U  =  7  m/s,  and  a  Bragg  wavelength  of  0.0  m,  corresponding  to  CgU  = 
0.069.  These  curves  are  roughly  at  6  =  ±  (C„/U)  cosy  for  the  starboard/port  wings,  so  that  0  varies  as  A  + and  U  - '. 
The  wing  angle  difference  Is  proportional  to  sln2vj.  For  purely  specular  reflection  (If  It  occurs),  one  can  see  only  one 
wing  or  the  other,  depending  on  The  specular  points  SP  where  the  Bragg  curves  cross  the  specular  curve,  the 
Bragg  scatterers  are  In  stationary  phase  with  the  ship;  for  those  special  conditions  one  can  expect  strcng 
signatures.  Doppler  offsets  are  Indicated  below  (SAR  traveling  from  right  to  left).  The  central  line  behind  the  ship 
represents  the  turbulent  wake  with  assumed  zero  Doppler  offset.  The  dashed  lines  give  the  positions  of  the  Bragg 
wings  without  Doppler  offset.  For  v  =  0  there  Is  no  Doppler  offset  for  the  ship,  and  a  slight  Doppler  offset  for  the 
wings  (to  the  right  for  the  starboard  wing,  to  the  left  for  the  port  wing),  as  Indicated  by  the  solid  V.  For  v/=90*  the 
Bragg  wings  coincide  with  the  turbulent  wake  (0  =  0)  and  are  not  offset,  but  the  ship  Is  offset  to  the  right.  At  the  SP 
point  of  v  =  8*,  the  starooard  Bragg  wing  has  the  same  Doppler  offset  as  the  ship,  whereas  at  <+»  =>  -8\  the  port 
Bragg  wing  has  the  same  Doppler  offset  as  the  ship. 


Figure  3-4.  The  wing  angles  6  as  a  function  of  the  look  angle  i  (positive  for  the 
starboard  wing,  negative  for  port). 


There  are  some  features  of  the  observed  V-wakes  which  are 


not  clearly  understood  in  terms  of,  or  perhaps  strain,  the  Kelvin 
wake  hypothesis.  We  list  these  below. 

(a)  Multiple  Crests 


SEASAT  Rev  407  shows  wakes  in  which  several  lines 
appear  to  be  visible  in  the  wake.  In  the  best  picture,  at  about  7  km 
behind  the  ship,  the  wavelength  between  the  two  outer  crests  is  about 
200  m,  and  between  the  two  inner  crests  about  100  m. 

The  version  of  the  Kelvin  wake  theory  described  above 
does  not  obviously  predict  multiple  crests.  Perhaps  they  result  from 
an  interference  pattern.  A  simple  two  point  source  model  will  give 
constructive  interference  at  sufficiently  small  wake  angles  if  the 
distance  between  the  two  point  sources  is  less  than  a  quarter  of  the 
transverse  Kelvin  wave  wavelength,  or  less  than  about  7m.  This  is  an 
exceedingly  small  number,  and  casts  doubt  on  this  idea. 

One  should  keep  in  mind  that  under  extreme  circumstances 
internal  wave  wakes  of  ships  can  be  seen.  The  very  strong  halocline 
in  Knight  Inlet  is  one  such  extreme  case.  Another  is  shallow  water; 
if  the  water  depth  is  only  a  few  times  the  mixed  layer  depth  very 
strong  internal  wave  wakes  can  be  produced. 
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V-wakes  are  quite  long,  up  to  perhaps  20  km 
lot  at  all  clear  that  30  cm  surface  waves 
survive  that  long,  even  in  conditions  of 
;ly  there  needs  to  be  more  study  of 
:t  wavelength  surface  waves,  in  various  wind 


A  20  km  long  wake  and  a  reasonable  guess  at  ship  speed 
suggests  perhaps  a  half  hour  persistence  time.  At  very  low  wind 
speeds  (less  than  2  m/sec  or  so)  the  lifetime  of  30  cm  surface  waves 
is  limited  by  viscosity.  The  e-folding  time  for  wave  amplitude  is 
t  =  (2  v  where  v  is  the  kinematic  viscosity  and  k  is  the 

surface  wake  wave  number.  For  30  cm  waves  this  time  is  about  15 
minutes.  Since  two  e-foldings  may  well  be  visible,  the  observed 
lifetimes  are  not  absurd,  provided  the  wind  speeds  are  really  low. 

The  lifetime  degrades  very  rapidly  as  the  wind  velocity  goes  up;  at  a 
wind  speed  of  6  m/sec,  for  example,  30  cm  waves  have  an  e-folding 
time  of  order  one  second  (see  Appendix  C.)  Therefore,  the  Kelvin 
hypothesis  is  only  tenable  provided  that  all  narrow  V-wake 
observations  occur  under  very  low  wind  speed  conditions. 

We  do  not  know  the  wind  in  pass  407 ;  the  dark  background  in 
the  image  suggests  it  could  be  quite  low.  In  Dabob  Bay,  wind  patches 
easily  visible  on  the  SAR  are  known  to  correspond  to  very  low 
winds.  Clearly  more  information  on  this  question  is  needed.  In  both 
Dabob  Bay  and  in  the  Gulf  of  Alaska,  the  narrow  V-wakes  disappear  in 
regions  of  wind,  and  no  wakes  are  visible  at  all  at  sea  states  of  3 


and  above. 


4.0  FUTURE  EXPERIMENTS 


We  believe  that  the  experiments  conducted  up  to  now  have 
demonstrated  that  the  narrow  V-wakes  seen  by  SARs  behind  surface 
ships  have  nothing  to  do  with  internal  wave  wakes.  We  also  believe 
that  they  have  nothing  to  do  with  the  classical  Kelvin  wake  of  a 
ship;  instead  we  have  suggested  that  the  observations  are  compatible 
with  a  simple  model  (described  in  the  previous  section  and  in 
Appendix  B)  in  which  the  near  axis  wake  is  ae  if  it  were  produced  by 
a  series  of  charges  dropped  behind  the  ship.  Exactly  what  this  model 
has  to  do  with  a  real  ship  remains  unclear;  presumably  the  near  axis 
wake  of  a  real  ship  is  governed  by  complicated  non-linear  and 
turbulent  flow  along  the  ship  sides  as  well  as  by  the  screw.  The 
"charges"  of  the  model  may  be  related  to  regions  of  turbulence  which 
are  sporadically  detached  from  the  hull.  To  understand  the  details 
of  all  of  this  will  be  very  difficult. 

Nevertheless  there  are  some  experiments  which  can  be 
performed  to  test  further  what  is  happening,  and  to  further 
distinguish  between  the  model  and  the  internal  wave  hyp/ "heels. 

These  lie  in  varying  Bragg  wavelength  of  the  radar.  In  the  internal 
wave  hypothesis,  the  radar  backseat  ter  is  caused  by  the  straining  of 
the  surface  due  to  an  internal  wave;  therefore  the  V-wake  angle  is 


determined  by  the  ratio  of  the  speed  of  the  ship  to  the  internal  wave 
speed.  The  V-wake  angle  should  therefore  be  independent  of  the  Bragg 
wavelength. 

In  the  "depth  charge  model",  on  the  other  hand,  the  V-wake 
angle  is  determined  by  the  ratio  of  the  group  speed  of  surface  waves 
at  the  Bragg  wave  number  to  the  ship  speed;  thus  changing  the  Bragg 
wavelength  will  change  the  V-wake  angle. 

The  Bragg  wavelength  can  be  changed  either  by  changing  the 
radar  frequency  fR  or  by  varying  the  radar  incidence  angle  i.  Since 
we  know  experimentally  that  the  narrow  V-wake  disappears  when  i  is 
too  large  (50°  or  more)  the  amount  of  variation  which  can  be  produced 
by  changing  i  is  limited.  Nevertheless,  changing  i  from  20°  to 
30°  will  change  the  group  speed  and  the  V-wake  angle  by  the  factor 
(sin  30°/ sin  20°)1/2  -  1.21  which  should  be  big  enough  to  measure. 
Larger,  and  therefore  more  easily  measured,  changes  can  be  made  by 
changing  fR.  However  it  is  important  not  to  use  frequencies  for 
which  the  corresponding  surface  waves  are  out  of  the  gravity  wave  and 
into  the  capillary  wave  regime,  since  the  physics  may  become  quite 
different  there. 

In  this  context  It  is  worth  noting  that  one  pair  of  ship 
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wake  images  taken  simultaneously  at  L-band  and  at  X-band  already 
exists. (*)  Both  wavelengths  show  narrow  V-wakes;  that  at  L-band  is 
much  longer  and  appears  wider  than  the  one  at  X-band.  The 


wavelength  ratio  is  7,  so  one  would  expect  the  ratio  of  wake  angles 

0 

to  be  vTT  The  observed  ratio  is  about  1  *  * _ 0  •  The  wake  at 

3.5  *  1 

L-band  is  much  longer  than  at  X-band.  Since  the  L-band  wake  extends 
off  the  picture,  the  wake  length  ratio  cannot  be  determined.  If  the 
lifetime  is  purely  due  to  viscous  dissipation,  the  ratio  should  be 
72.  The  length  of  visible  track  at  X-band  seems  to  be  of  order  a 
kilometer.  The  ship  speed  can  be  determined  from  a  clear  transverse 
wave  modulation  of  the  wake;  it  is  about  18  knots.  Thus  the  X-band 
surface  waves  (whose  wavelength  is  4  cm)  must  therefore  have  lived 
some  100  seconds.  This  is  an  uncomfortably  long  time,  and  reinforces 
our  earlier  remarks  about  wake  lengths. 


(*)  These  were  taken  in  July  1978  using  an  aircraft  SAR.  See  R. 
Hammond,  et.  al.,  NOSC  interim  report  on  "Joint  U.S.  -  Canadian 
Ocean  Wave  Investigation  Report,"  p.44. 


APPENDIX  A 


SPECKLE  NOISE 

Since  detecting  weak  hydrodynamic  modulation  of  SAR  return  is 
a  central  issue  it  is  worth  examining  the  intrinsic  SAR  noise,  the 
dominant  one  being  speckle  noise. 

The  change  of  radar  cross-section  produced  by  particular 

hydrodynamic  processes  can  not  be  predicted  quantitatively  nor  is 

absolute  calibration  of  SARs  straightforward.  Consequently,  it  is 

convenient  to  reference  modulation  of  return  to  the  average  background 

return.  Thus  if  I  is  a  SAR  intensity,  <1>  is  an  average  which 

suppresses  SAR  noise  but  passes  any  modulation  <6I>  of  the  intensity 

and  <I>0  is  ^  average  <I>  of  the  background,  then  the  most 

convenient  measure  of  signal  strength  is  the  modulation 

index  <6I>/<I>  .  Unfortunately,  receiver  noise  makes  it  difficult  to 
o 

establish  zero  intensity  so  that  measured  intensity  is  often  in  error 

by  a  constant.  Consequently,  reported  intensity  ratios,  particularly 
at  low  signal  levels,  may  be  in  error. 


A-l 


If  the  SAR  return  is  dominated  by  first  order  Bragg 
scattering  the  intensity,  I,  of  a  single  radar  look  should  be 
distributed  as  a  chi-squared  variable  with  two  degrees  of  freedom;  in 
particular  the  variance  of  I  would  equal  the  squared  mean 

'2  2 
<1  >  =  <I>  > 

Any  hydrodynamic  modulation  must  be  strong  enough  to  stand  out  above 
this  speckle  noise. 

In  the  Seasat  SAR  images  digitally  processed  by  JPL,  four 
radar  views  separated  along  the  satellite  track  are  averaged 
incoherently  to  produce  an  approximately  12m  by  12m  pixel.  JPL  images 
are  typically  digitized  In  units  of  amplitude.  A,  which  is  the  square 
root  of  the  energy  unit  I.  If  the  energy  in  N  independent  looks  are 
averaged  the  resulting  intensity  and  amplitude  obey  the  relations 

<l"2>  -  |  <I>2  (A-l) 

<k'Z>  -  -  N  ]  <A>2  (A-2) 

r  (N+i/2) 

We  have  examined  four  100-pixel  regions  from  a  JPL  version  of 


the  Seasat  Rev  407  image  near  Cape  Canaveral  in  which  a  wake  with 


banded  structure  has  been  identified.  The  four  regions  were  examples 
of  unmodulated  background  found  near  this  wake.  Analysis  indicates 
little  interpixel  correlation  except  that  neighbor  pixels  separated  in 
the  along-track  direction  had  correlation  0.6  ^  0.1.  Robust  intensity 
and  amplitude  variance  values  were  found: 


^  >  TT.W.7T 

/A '2s  _  1 

^  '  m+iT 


<i>2 

<A>2 


(A-3) 

(A-4) 


Since  the  four-look  pixels  are  correlated  in  the  azimuth 
direction,  so  would  the  four  individual  looks  in  an  average  be 
correlated.  Since  a  pixel  is  an  average  over  less  than  four 
independent  looks  the  speckle  noise  is  more  energetic  than  if  four 
independent  looks  had  been  averaged.  In  fact  (A-3,  A-4)  are 
consistent  with  (A-l,  A-2)  with  N  =  3.  The  observed  azimuth 
correlation  is  somewhat  larger  than  is  easily  explained  by  this  slight 
reduction  in  the  effective  degrees  of  freedom. 


In  summary,  the  speckle  noise  in  JPL  processed  Seasat  SAR 
images  is  reasonably  consistent  with  a  first  order  Bragg  scattering 
mechanism  and  a  modest  correlation  of  the  fundamental  looks  from  which 
the  four-look  pixels  are  obtained. 


APPENDIX  B 


FIRST  ORDER  BRAGG  CROSS  SECTION:  LIMITING  CASES 


Suppose  we  model  the  surface  ship  as  a  known  but  arbitrary 


pressure  distribution  pQ(x,  t)  on  the  water  surface.  Then  the 


(linearized)  surface  displacement  s(x,  t)  produced  by  this  pressure 
is 


s(x,t)  -  / 


d2k 


i£«x 


(2tt) 
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/  dt '  j  d2x 


/  k  -ik  •  x  +,+'  . .  exp  < 
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The  first  order  Bragg  backscatter  cross  section  is 


,2  H.  -2ikH  '  x 


proportional  to  o(t)  ■  2ikv  /  d  x  e  “  s(x,t)  where 

the  ky  and  ^  are  the  vertical  and  horizontal  radar  wave  numbers, 
and  the  integral  is  taken  over  the  illuminated  area.  For  simplicity 
let  us  suppose  the  illuminated  spot  is  defined  by  a  Gaussian  cutoff 


.2 

-  a  (x-x  )  + 

e  centered  on  a  point  x  with  a  spot  size  of  dimension 
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1/a  . 
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Then  we  find 
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Here  C 
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P*- 


waves  at  the  Bragg  wave  number 


is  the  group  velocity  of  surface 

2V 


This  expression  is  valid  for  an  arbitrary  pressure 
distribution  p^x.t).  Let  us  now  look  at  some  limiting  special 


cases. 


A  moving  point  source. 


Here  pQ(x,t)  =  pQ  -  tjt)  where  Tj  is  the  ship 


velocity.  Then 
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Aside  from  an  overall  phase  exp  (-2i  ^  •  tft),  o  is  a  function  only 


■j  V 


An  (approximate)  evaluation  of  (B— 3)  yields 
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Evidently  o  will  be  exponentially  small  unless  both  of  the 
conditions 


and 
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(B-5b) 


(B-5c) 

(B-5a) 


are  satisfied.  When  these  are  satisfied,  we  have 


TT 

/  a 


l 

U 


(B-6) 


The  pair  of  conditions  (B-5a,  b)  or  (B-5c,  d)  indicate  that 
the  radar  sees  only  waves  having  the  Bragg  wavelength  and  with  crests 
oriented  perpendicular  to  the  radar  look  direction. 


Let  us  choose  the  ship  to  be  moving  in  the  x  direction. 
Then  condition  (B-5a)  says  that  a  return  exists  only  at  the  wake 
angle  3  satisfying 


tan  0  = 


C  C  cos  t|> 

gy  =  -  g 

U  -  C  U  -  C  sin  ii> 

gx  g 

where  ^  is  the  angle  between  the  satellite  track  and  the  ship  track; 

or,  for  C  «  U  ,  and  therefore  0  small, 

g 

C 

q  M  C08  ^  > 

(the  other  arm  of  the  V-wake  comes  from  (B-5c))  which  is  the  usual 
Bragg  selected  wake  angle.  Condition  (B-5b),  however,  also  requires 
that 

2 V '  ^ 


This  is  the  Kelvin  wake  condition;  it  requires  that  the  radar  be 
oriented  perpendicular  to  the  Kelvin  wake  wave  crests.  This  will  be 
the  case  if 


gx 


> 


therefore  only  for  a  single  radar  aspect  angle  relative  to 
wave  is  there  a  return.  (See  Figure  B-l.)  In  particular, 
only  one  arm  of  the  Kelvin  wake  will  be  visible. 


the  wake 
therefore, 


B. 


A  line  of  fixed  point  sources 


Next  let  us  take  the  source  pressure  to  be 

p  v,x ,  t)  =  p  Z  S“(x  •  3  n  t)  6(t  -  nx)e  n 
ro  ro  n 

representing  a  line  of  stationary  (in  the  water)  instantaneous 
incoherent  sources.  (This  is  the  simple  depth  charge  model  described 
in  Section  III).  Now  we  find 

2i  kv  p  /2k~  -2i£  •  £nx  i<|> 

a  = - -  /  -SL  z  e  H  e  n  0  (t-nx) 

P  /  g  n 

,  -i/2gk  (t-nx)  -a  [x  -  tfnx  -  £  (t-nx)]2 

If  H  '•o  g 

Ti  'e  e 

-  (Je  *  -*8*  *  /<*?)  <“-7) 


Suppose  the  spacing  between  sources  is  big  enough  so  that  only  a 


single  term,  say  n  *»  no,  in  the  sum  contributes  for  a  given  spot 

location  x  and  time  t.  Then 
o 


/  8 


P  /  g 

and  we  must  have  the  constraint  x  ±  £  t 

0  2 


(ft  ±  £  )n  t  « 
g  o 


(*)  Keep  in  mind  that  here  has  different  dimensions  from  the  PQ 
in  the  limit  described  in  Eq.  B-6. 
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Thus  we  obtain  returns  at  wake  angles 


tan  0  *  ?  — - - 

U  ±  C  sin  i|i 

gx 

or,  for  small  angles, 

C 

8  *«  ±  cos 

The  wake  angles  allowed  are  the  same  as  in  Case  A  (constraint  B-5a  or 
B-5c)  but  there  is  now  no  constraint  on  the  radar  orientation 
relative  to  the  wake  direction;  (see  Figure  B-2)  there  is  a  return 
for  any  direction  of  ^  . 

If  the  spacing  between  sources  is  much  less  than  the  spot 
size,  and  if  the  phases  <t>n  are  all  the  same  so  that  the  sources  are 
coherent,  then  this  situation  reduces  to  Case  A,  the  normal  Kelvin 
Wake.  Any  intemediate  case,  or  a  case  with  any  other  kind  of 
pressure  distribution  in  space  and  time,  is  described  by  Eq.  (B-2). 

Under  some  circumstances,  second  and  higher  order  Bragg 
scattering  can  be  Important.  The  second  order  Bragg  cross  section  is 
calculated  in  other  JASON  reports. 
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APPENDIX  C 


NOTE  CONCERNING  THE  PERSISTENCE  OF 
A  PATTERN  OF  SURFACE  GRAVITY  WAVES 


I .  INTRODUCTION 


The  imaging  of  surface  ship  Kelvin  wakes  by  a  SAR  raises  an 
interesting  question:  How  long  will  a  given  "pattern"  of  surface 
waves  persist?  Under  conditions  of  light  wind  and  low  sea  state  the 
Kelvin  wake  of  a  vessel  may  be  seen  tc  perist  astern  for  some 
kilometers.  Under  more  robust  sea  conditions  the  wake  persistence 
seems  to  be  greatly  reduced.  In  this  note  we  shall  discuss  aspects 
of  the  decay  of  such  a  "pattern".  This  is  a  complex  phenomenon  and 
cannot  be  exhaustively  addressed  here. 

Deterministic  mechanisms,  such  as  linear  wave  propagation 
and  the  interaction  with  known  large  scale  currents,  can  distort  a 
wave  pattern.  This,  however,  is  in  principle  predictible  and  will 
not  be  considered  as  "pattern  decay".  Stochastic  mechanisms,  such  as 
interaction  of  the  pattern  with  wind  and  ambient  sea,  do  lead  to  a 
genuine  decay  of  the  wave  pattern. 


We  shall  be  concerned  only  with  patterns  formed  by  gravi  ty 
waves.  Capillary  wave  decay  is  expected  to  occur  at  time  scales  too 
short  to  be  of  interest  here. 

There  are  evidently  several  mechanisms  that  contribute  to 
wave  "pattern  decay".  Viscous  dissipation  seems  to  be  the  simplest 
of  these.  Wind-wave  interaction,  while  physically  complex,  is 
phenomenologically  modelled  as  a  simple  exponential  pattern  decay. 
These  mechanisms  are  reviewed  in  Section  3. 

Nonlinear  hydrodynamic  interactions  of  the  "pattern"  with 
ambient  surface  waves  provide  other  decay  mechanisms.  In  the 
terminology  of  wave-wave  weak  interaction  theory,  both  "three-wave" 
and  "four-wave"  interactions  must  be  considered.  It  is  well-known 
that  the  three-wave  interactions  do  not  admit  resonant  frequency 
conditions.  We  shall  see  that  these  lead  to  only  partial  pattern 
decay.  Only  if  this  decay  brings  the  signal-to-noise  level  below  ttfie 
threshold  for  detectibility ,  can  we  consider  the  pattern  destroyed. 
The  four-wave  Interactions  act  on  a  slower  time  scale  than  do  the 
three-wave  interactions,  but  lead  to  total  decay  of  the  pattern 
(according  to  weak  interaction  theories). 

The  role  of  three  .ind  four-wave  interactiosn  is  discussed 
in  Sections  4,  5,  and  6.  In  Section  4  we  discuss  the  decay  with 
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respect  to  an  amplitude  matched  filter  in  the  space  time  domain.  In 
Section  5  we  describe  the  evolution  of  the  power  spectrum  of  the 
pattern  due  to  three-wave  interactions.  In  Section  6  we  make  some 
estimates  of  the  spectral  evolution  resulting  from  four -wave 
interactions. 

2.  DESCRIPTION  OF  THE  DECAY  PHENOMENA 


In  a  specified  "rectangular  area  of  ocean"  we  use  a  Fourier 
expansion  for  the  surface  wave  displacement  C(x,t)  (  x  «*  x,y  ,  a 
vector  in  the  plane  of  the  quiescent  ocean  surface): 

5(x,t)  «  -Im[  Z(x,t)] 

z(x,t)  -  exp[i(£*x  -  w£t)]  .  (C-l) 

fst 

Here  m  =  (g£)^  .  In  the  approximation  that  the  waves  are  linear 

Xt 

and  not  forced,  the  b.  are  constants.  In  general,  b _ ( t )  is  time 

X*  Jc 

dependent . 


We  shall  suppose  that  (C-l)  refers  to  the  ambient  sea.  For 
the  "pattern  waves"  we  take  as  a  special  case  of  (C-l) 

C  (x,t)  ■  -Im[z  1 

p~  1  pJ 
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(C-2) 


Zp  "  kE  \  exP[i(£  *  2  “  \fc)  1 

*■  *+ 

We  shall  suppose  that  the  are  amplitudes  (they 

determine  the  pattern)  at  some  reference  time,  say  t  *  0.  The 
ambient  sea  will  here  be  presumed  as  linear,  so  the  b  are 

Ar 

constants.  Over  an  ensemble  of  realizations  of  the  ocean  surface  the 
b.  will  be  considered  to  be  uncorrelated  Gaussian  variables.  The 

Xf 

/v 

ensemble  averaged  spectrum  ¥  of  vertical  displacement  is 

/  K£)<i2a  -  j'<  Ib^l2  >  (C-3) 

In  this  paper  we  shall  use  a  Pierson-Moskowitz  (1964)  spectrum  for 
the  ambient  sea: 

»<*>  -  (-J)  «p[-  C%)  ]  G(6)  , 

i 

n  «  4  x  10"3  ,  Zq  -  8.4  /  w2  (04) 

Here  W  is  the  wind  velocity  expressed  in  m/sec,  Zq  is  in  m  *  , 
and  3  is  the  angle  between  Z  and  W  .  We  shall  specify  models  of 
G(9)  later,  but  now  note  the  normalization 

tt 

/  G(9)  d0  ■  1  . 

-it 


04 


(05) 


For  the  pattern  waves  we  introduce  the  Wigner  (1932) 


correlation  function 


r 

p 


(£>5.0 


7  <  zp(x  -  Z*  (I  +  £/2.t)  > 


(C-6) 


Here  <•••>  represents  an  ensemble  average  over  the  ambient  sea,  as 

* 

in  (C-3).  The  spectrum  of  vertical  displacement  is 


¥  (k,x,t)  * 


ik*r 

A.  fS. 

e 


d2r/ (2ir)2 


(C-7) 


The  spectrum  of  wave  action  F^  and  the  energy  spectrum  (both  per 
unit  area  of  ocean)  are  related  to  (G— 7)  by  the  expressions 


(C-8) 


Here  pQ  is  the  density  of  sea  water. 


form 


The  Hasselmann  (1967)  evolution  equation  for 


F  is  of  the 
P 


3F 

-5~£+A«VF  -  s.  +  s  +  s 

dt  xp  ni  v  w 


(C-9) 


It  is  assumed,  of  course,  that  varies  slowly  over  distances 
comparable  to  pattern  wavelengths  of  interest. 
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Here 


j  ■  W  * 

and  S  S  ,  and  S  represent  the  effects  of  wave-wave 
n£*  v*  w  r 

interactions,  viscosity ,  and  wind,  repectively.  In  (2.9)  we  are 
omitting  interaction  with  large  scale  surface  currents. 


3.  EFFECTS  OF  VISCOSITY  AND  WIND 


Viscosity  will  lead  to  exponential  damping  of  the  pattern 
waves.  In  the  notation  of  (C-9)  this  is  expressed  as  (Phillips, 
1977) 


S 

v 


-0  F 
v  p 


9 


3 


v 


4vk2 


> 


(C-10) 


<%  -6  2 

where  v  —  1.1  x  10  m  /sec  is  the  kinematic  viscosity  of  sea 
water.  The  decay  time 


v¥) 5  C 


(C-li) 
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is  shown  as  a  function  of  wavelength 


X 


in  Figure  C-l. 


The  rate  of  wind  induced  growth  for  surface  waves  has 
recently  been  reviewed  by  Plant  (1982).  He  concludes  that  for  the 
growth  rate  of  small  amplitude  waves, 

S  -  -8  F  (c-i: 

w  w  p 


in  (C-9).  Here 


w 


2  2 
0.04  u^k 
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6/\ 


(c-i: 


where  u*  is  the  friction  velocity  of  the  wind  and  0  is  the  angle 
between  k  and  w  . 

#s»  ^ 


It  seems  plausible  for  our  purposes  to  interpret 


To  <">  5  £ 


as  a  pattern  decay  time. 


(C-l 


To  evaluate  ^.C-14)  we  use  the  analysis  of  Garratt  (1977)  to 
relate  u*  to  the  wind  speed  W.  The  resulting  time  T^(w)  shown 
in  Figure  (C— 2)  as  a  function  of  W  for  several  wavelengths  X  and 
angle  0=0. 


We  note  from  Figures  (C-l)  and  (C-2)  that  for  W  >  2  m/sec 
that  Tj(w)  <  Td(v)  when  X  >  o.l  m  • 

4.  DECAY  FROM  A  PREDICTED  PATTERN 


One  simple  descriptor  of  pattern  decay  is  given  by  the 
correlation  function 

r  (t)  =  <  P{;  (x,0;t}  C  (x,t)  >  CC— 15) 

r  r 

Here  t  =  0  is  considered  to  be  the  reference  time  at  which  the 
pattern  is  first  observed.  The  quantity 

P  {  Cp(x,0);t)  | 

represents  the  predicted  field  at  t,  given  the  initial  field 
4p(x,0)  .  The  prediction  Is  to  be  obtained  using  known 
deterministic  phenomena  which  may  distort  the  pattern.  For  example, 
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(soas)  Pj. 


if  linear  wave  propagation  is  the  only  deterministic  phenomenon,  then 


P  {  ?  (x,0);t)  }  =  -Im{kI  Bk(0)  exp[i(k*x-o^t)  ]  }  (C-16) 

*  ^  M 


where  Bk(0) 


is  the  value  of  ^  at  t  =  0. 


In  this  paper  we  shall  assume  that  (C-15)  applies  in 
evaluating  (C-15).  Then  it  is  sufficient  to  calculate  the  set  of 
correlations 


Ykk,(t)  =  <  Bk(t)  Bk/(0)  >  -  <  Bk(t)  >  Bk,(0)  (C-17) 

(\«V  A/  (V  /V  ^ 

We  shall  refer  to  a  given  Fourier  amplitude  in  the  pattern  as  a  "test 
wave" . 


Under  certain  conditions  we  might  expect  <  Bk  >  to  be 

A/ 

determined  from  a  Langevin  equation  of  the  form 


ft  <  Bk  >  -  -V(k)  <  Bk>  ,  (C-18) 


where  v(k)  is  the  Langevin  "rate  constant".  In  this  case  we  would 
evaluate  (C-17)  as 


Ykk,“  »k(0)  B  .  ,(0)  exp[-v(k)t]  . 


(C— 19) 


We  shall  see  that  four-wave  interaction,  with  frequency  resonance, 
lead  asymptotically  to  equations  of  the  form  (C-18)  and  (C-19). 

As  a  first  illustration  of  (C-15)  and  (C— 16)  we  consider  the 
advection  of  a  test  wave  "k"  due  to  the  surface  current  U(x,t)  of 
the  ambient  waves.  This  current  ie  of  the  form 

JJ(x,t)  =  E(i)  l  w  (b  exp[i(fc‘x-m  t)]  -  c.c.}  (C-20) 

ft  ^  ft  * 

using  the  notation  of  (C-l).  We  have  assumed  that  the  ambient  waves 

are  linear,  so  the  b  are  constant.  We  shall  also  assume  that 

a 

£  <  <  k  ,  (C-21) 


or  that  the  wavelength  of  the  test  wave  Is  small  compared  with  that 
of  those  ambient  waves  which  are  of  most  importance  for  pattern 
decorrelation.  This  assumption  can  be  tested  for  verification  when 
we  evaluate  the  decay  using  the  spectrum  (C-4). 

If  the  test  wave  is  advected  with  the  local  velocity  U  we 
use  (C-21)  to  write 


t 

\(t)  -  Bk(0)  enpl-ik  •  /  U(t'  )dt'  )  C  C— 22 ) 

O 


c-u 


Using  the  assumed  Gaussian  distribution  of  the  b^'s  and  the 
relation  (C-3)  we  obtain 

<  B  (t)  >  =  Bfc(0)  exp[-D3(k,t)]  , 

<v 

where 

D3(k,t)  -  /  d2i  (k»i)2  [l  -  cos(a>At)]  ¥(i)  . 

For  small  t  this  becomes 

D,(k,t)  -  o2t2/2  ,  t  <  <  1  , 

3  ~  £o 

where 

o2  “  <  (k*Ji)2  >  -  I  d2*(k«l)2  u^C*)  . 

For  t  +  •  ,  on  the  other  hand,  we  have 

D3(fc>«)  -  /  d 2l  (k/I)2*(A)  . 

We  note  that  for  the  spectrum  (C-4)  the  assumed  condition  (C— 2 1 
appears  to  be  valid  for  the  evaluation  of  D-  . 


When  (C-25)  is  valid  we  may  define  an  e-folding  decay  time 


as 


Td(3)  =  /2  /  o 


(C-28) 


More  generally,  we  may  define  the  decay  time  with  the  relation 


D3(k,Td(3)  )  -  1  .  (C-29) 

This  equation  may,  or  may  not,  have  a  solution.  When  there  is  no 
solution,  the  advection  mechanism  is  ineffective  in  destroying  the 
"pattern".  The  single  e-folding  condition  (C-29)  is  arbitrary  and 
one  may  wish  to  define  the  decay  time  to  correspond  to  several 
e-foldings. 

A  more  systematic  and  formal  description  of  the  decay  of  (C- 
17)  can  be  obtained  using  the  method  of  Van  Kampen  (197A).  To 
develop  this  we  begin  with  the  formulation  of  Watson  and  West 
(1975).  Their  equation  of  motion,  when  linearized  in  the  pattern 
amplitudes  B^  ,  is  of  the  form 


\  -  T2(B,b)  +  X3(b,b)  +  .  .  . 


(C— 30) 


Here  T,?  and  T3 


are  quadratic  and  cubic,  respectively  in  wave 


C — 1 4 


amplitudes.  If  we  ignore  all  but  the  term,  (C-30)  takes  the 
form 

\-pEAk\  •  <c-31> 

where 

E  k 

\  '  V-J-fc  l2rE*  \  exPt1(\-“p+“l)t] 

k-A  * 

+  Tp  ~  ~b_^  exp  [ i ( \~Wp~a)^) t  ]  }  (0-32) 

The  coefficients  T  are  defined  in  Appendix  A  of  Watson  and  West 
(1975). 

The  Van  Kampen  equation  is 

a?  <  V  ■  -  Wt>  <  \  >  •  (c'33) 

where 

t  £  k 

K(t)  -  -  /  dr  <  Ak  (t)  A  ~(t  -  t)  >  .  (C-34) 

Q  IL 
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Evaluation  of  this  is  straight  forward,  assuming  that  the  b  are 

X 

Gaussian  and  the  l  <  <  k,  p  .  The  result  is 

K(t)  =*  /  d 2 l  (k*i)2  sin(wAt)f(4)  (C-35) 

Integration  of  (C-33)  leads  to  the  expressions  (C— 23)  and  (C-24). 

Yue  three-wave  interaction  model  (C-31)  is  the  equivalent  to 
the  simple  udvection  model  described  by  (C— 22) •  The  fact  that  our 
decay  function  does  not  increase  indefinitely  with  time  reflects 

the  lack  of  three-wave  resonance.  To  find  true  decay  to  a  vanishing 
pattern  amplitude,  we  include  the  effects  of  7^  in  (C-30).  This 
contains  four-wave  interactions.  Higher  order  terms  arising  from 
T^  also  contribute  four-wave  interactions  and  must  also  be 
Included. 


Taking  account  of  the  four-wave  Interactions,  we  again  have 
an  equation  of  the  form  (C-31),  but  with  the  definition 


kn 

*v  ^  r  ** 


X 


exp  [i(<Y» 


01  +  0) 

P  a 


(C-36) 


C-16 


The  coefficients  C  here  are  defined  in  Appendix  B  of  Watson  and  West 
(1975)  [see  specifically  Eq.  (47)  of  this  reference].  Some  rapidly 
oscillating  terms  have  been  dropped  from  (C-36).  These  terms  do  not 
contain  frequency  resonance  and  do  not  lead  to  a  true  pattern  decay. 

Evaluation  of  (C— 34)  using  (C-36)  now  gives  the  decay  rate 


v(t)  -  Re  [K(t)] 

■  m  !  -l2MZA*U>«8>t“nfe  •  U+  n)  ]2 

v  atn  (90 

X  3 


where 


3  - 


-  0)  + 
n 


*  **  **  #v» 


■  V  wn  +  £(y  •  («  -  4) 


(C-38) 


and 


c(k)  . 

As  t  +  “»  ,  we  obtain 
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v(«)  “  /  d2jl  d2°  Vn[^  *  +  E>  l2  5<S)  (C- 


We  see  that  asymptotically  the  four-wave  resonance  leads  to  a 
Langevin  decay  law  of  the  form  (C-19). 

For  finite  times,  we  have  the  relation 

<  Bk(t)  >  -  Bk(0)  exp[-D4(k,t)  ]  (C- 

A/ 

where 

t 

D4(kt)  »  /  v(t')dt' 
o 

"  ils  / dZw2n  “a  k*(&  +  n)  l2 

X  [l  -  cos(Bt)  ]  /  e2  (O 

As  t  +  “  ,  we  obtain 

B4(k,t)  -  v(*)t  .  (C- 

To  describe  the  implication  of  (025)  and  (041)  we  first 


assume  an  isotropic  ambient  spectrum: 


G(  8)  =  1  /  (2ir) 


(C-43 


Although  not  realistic  as  a  wind  wave  spectrum,  (C-43)  leads  to 
reasonably  accurate  decay  rate  because  of  the  integration  over  wave 
angles.  A  more  realistic  spectrum  will  be  considered  later. 


Using  (C-43)  we  can  readily  express  (C-25)  in  form 


D3(k,t)  =  (  £  )  (  k  /  £q]2  J3(s) 

00 

J.(s)  =  /  [l  -  cos(/x  s)  ]  dx  /  X3  , 
J  1 


(C-44; 


where 


8  fiV 

o 


(C-45! 


and  we  have  replaced  the  exponential  factor  in  (C-4)  by  a  simple  cut 
off  at  l  »  .  Similarly,  we  may  evaluate  (C— 41)  as 


“  (  %  )  (k  I  A0)  » 


In  Figure  C-3  we  show  the  limiting  decay  (C-27)  for 

three-wave  interactions  for  several  wavelengths  X  =  2n  /  k  as 

functions  of  wind  speed  W.  When  D  (k,®)  is  small  enough  that  the 

3  ~ 

pattern  cannot  be  considered  destroyed,  the  three-wave  mechanism  is 
ineffective.  We  have  arbitrarily  chosen  here 

D3(k,«)  >  1  (C-4 

as  the  condition  of  pattern  destruction  by  three-wave  interactions. 

We  now  define  a  decay  time  1  as 

d 

T,  ■  smaller  of  T.,  T.  , 
a  3  4 

D3(k,  Tj)  -  1  , 

D4(k,  T4>  -  1  .  (C-4 

In  Figure  C-4  we  show  as  a  function  of  wind  speed  for  several 
wavelengths  X  .  The  dashed  lines  represent  an  interpolation  made 
from  T.  (upper  region)  to  t  (lower  region). 


I 


cy.' 

> 


)  due  to  three-wave  interactions  is  shown  as 
ral  wavelengths. 


/  / 


A  more  realistic  spectrum  than  CC— 43)  is  that  of  Tyler  et 
al.  (1974): 

F( 9)  -  cosS(6)  /  L(S)  , 

L(S)  =  2/?r(jS+-)/r(|-S  +  l)  ,  (C-49) 

where  S  is  a  function  of  k. 

Using  (C-49)  we  evaluate  (C-26)  as 

oo  o  o 

o  =  ng(kz  /  [0»8  cos  ct  +  0.2  sin  a  ]  , 


a  *  0.4  (W  /  X)  (cos^a  +  0.25  sin^a)^  .  (C-50) 

Here  a  is  the  angle  of  k  with  respect  to  W  .  Use  of  (C-43),  on 
the  other  hand,  would  replace  (C-50)  by 

o  =  0.3  (W  /  X)  . 

The  consequence  of  using  the  spectrum  (C-49)  does  not  seem  very 
significant. 
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superposition  of  plane  waves,  can  readily  be  analyzed  by  the  method 
presented  in  this  section. 


The  wind  velocity  vector  W  =  iW  is  assumed  to  be  directed 
parallel  to  the  x-axis.  The  b  's  are  supposed  to  have  a  Gaussian 
distribution  and  the  ambient  spectrum  is  assumed  to  be  described  by 
(C-4).  The  expression  (C-6)  can  be  evaluated  analytically.  It  is 
consistent  with  the  condition  (C-21)  to  assume  that  ti  <  <  1  ,  so 

we  obtain 


— i]3  *r 

rp(r,t]  -i  |B£(0)|2  A  ~  exp{-  /  d2*  ¥(A) 

x  (*  *  *)2(a  *  rfC1 '  cos(wA0  ]  }  ♦  <c-5'+ 

Equation  (G-7)  may  also  be  evaluated  analytically  using  the 
expression  (C-54).  The  result  gives  the  pattern  spectrum  as  a 
function  of  time: 

fpfet)  =  |b£(0)  I2  /  [8.  ] 

x  exp[  -  (K^cosij)  -  k  sin /  (4c.) 

•  x  y  i 

-  (Kx  sin<j>  +  Ky  cos <t>)2  /  (4c2)  ]  (C-55) 


C1 

=  cos  <j> 
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2 , 
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cosi|)  b 

tan  2 $ 
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r  ’ 
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’  *)2[‘ 
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a2 
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•  *)2(n  • 

■  i)2[> 

-  cos(wAt)  ] 

b 

=  /  d 2l 

*(»)  (i  •  ii)(i  • 

*)(e  • 

i)[l  -  cosCm^t 

(C— 56) 

It  is  straightforward  to  evaluate  these  expressions  for  the 
spectral  form  (C-49).  This  is  tedious  and  seems  overly  elaborate  for 
our  purpose.  We  therefore  consider  only  the  "peaked"  spectrum 
corresponding  to 


G(0)  -  6(9) 


(C-57) 


and  the  isotropic  spectrum  (C— 43) 


For  the  peaked  spectrum  (G— 57 )  we  may  simplify  (C-55)  to  the 


form 


|b£(0)|2 

yk>t)  -  — — —  pj 

4/  ifa,  y  y 


x  exp[-  K2  /  (Aa^  ] 


(C-58) 


In  this  case 


a 


1 


(C-59) 


where 


h  -  n  /  y  exp[-  {lQ  /  i)2]  [l  -  coB(w£t)]  ( C— 60) 


For  the  isotropic  spectrum  (C— 23 )  we  suffer  no  loss  of 


generality  in  setting  p  ■  0  . 


I  MO)  | 2 

yfe.t)  — — - 


Then 


exp  [-  (kx-  px)2  /  (4ax)  ] 


X  exp[-  k*  /  (4a2)  ]  . 


(C— 61 ) 


Now, 


a!  -  |  £  h  •  V  £  4  h  •  ”  ■  0  <C-62) 

where  h  is  given  by  (C-60). 

The  integral  (C-60)  requires  a  cut-off  at  large  i  and 
depends  logarithimically  on  this  cut-off.  Thus,  the  scale  separation 
condition  (C— 2 1 )  is  not  automatically  satisfied  by  the 
spectrum  ( C— 4 ) ,  but  must  be  imposed.  For  numerical  evaluation  we 
have  done  this  be  inserting  a  factor 

[i  +  (2i  /  px)z  rl 

into  the  integrand  in  (C-60).  In  spite  of  his  arbitrary  imposition 
of  a  scale  separation,  our  spectral  evolution  equations  appear  to  be 
of  some  interest. 

As  was  the  case  with  the  decay,  as  described  by  (C-27),  the 
spectral  spreading  reaches  an  asymptotic  limit  for  large  t.  To 
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illustrate  our  model,  we  consider  the  peaked  model  (C-58)  for  the 
case  that  p^  =  0  .  The  spectrum  is  shown  in  Figure  C-5  for 
X  =  2ir  /  p  =  0.1  m  and  w  =  1  and  10  m/sec.  In  Figure  C-6  we  shown 
the  spectrum  for  X  =  0.4  m  and  w  =  10  and  20  m/sec.  The  spectrum 
described  by  (C-61)  is  similar,  but  evolves  in  both  horizontal 
dimensions . 

For  these  examples  the  asymptotic  variance  in  the  wavenumber 
is  much  less  than  the  initial  wavenumber  of  the  wave  train. 
Nevertheless,  this  spreading  can  decorrelate  the  wave  phase  coherence 
for  a  wavetrain  of  many  wavelengths. 


6.  DIFFUSION  IN  WAVENUMBER  SPACE 


The  four-wave  interactions  lead  to  an  irreversible 
re-distvibution  of  wave  energy  in  wavenumber  space.  This  is 
described  by  the  term  Sn^  in  (C-9).  In  the  general  case  this  is  a 
very  complex  phenomena  and  beyond  our  present  scope  to  discuss.  For 
the  case  that  the  scale  separation  condition  (C-21)  can  be  assumed 
Sn£  reduces  to  a  simple  diffusion  mechanism 


a  F. 


ij  3k. 


(C— 63) 
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Weak  interaction  theory  leads  to  the  expression 


D 


ij 


Z  3U  Z  dU 

=  i  <  f  k  •  ~  dt '  f  k  »  ~dt"  > 
t  *  ~  3x .  1  ~  3x . 


j 


t  00 

(C-64) 


Here  U  is  given  by  (C-20).  For  our  application  with  the  spectrum 

(C-4),  the  scale  separation  condition  must  be  imposed  by  introducing 

a  cut-off  at  wavenumber  greater  than  some  value,  say  l  .  For 

may 

this  reason,  only  a  brief  qualitative  discussion  of  (C-63)  and  (C-64) 
seems  justified.  A  scalar  diffusion  coefficient, 


appears  sufficient  then,  to  characterize  '<  ,ie  diffusion.  For  the 
spectrum  (C-49)  we  obtain 

D  =  3  x  10  ^  k^  w  ,  0)  =»  /~g  i  .  (C-65) 

c  c  ma 

y 

The  variance  Ak  in  an  initial  wave  of  wavenumber  k  is  estimated  at 
time  t  as 


A  k  »  /~"bt 


-  2  x  10 


-2 


/"(*)  t 
c 


(C-66) 
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If  we  assume  ~  1  sec  *  as  reasonable,  (C-66)  implies 
times  of  the  order  of  10^  seconds  for  significant  pattern 
distortion.  This  is  not  incompatible  with  the  four-wave  time  scales 
of  Figure  C-4. 

7 .  CONCLUSIONS 

We  have  seen  that  the  decay  of  a  wave  pattern  on  the  ocean 
surface  is  a  complex  phenomena.  A  number  of  mechanisms  may 
contribute  to  this  with  many  orders  of  magnitude  variation  in 
lifetimes,  depending  on  pattern  wavelengths  and  environmental 
conditions . 
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APPENDIX  D 


CORRECTIONS  TO  SAR  ESTIMATES  OF  V-WAKE  ANGLES 


1.0  Introduction 


Interpretation  of  synthetic  aperture  radar  (SAR)  images  of 
ship  wakes  requires  a  detailed  knowledge  of  some  aspects  of  wake 
image  shape.  For  example,  the  included  angle  of  the  wake  allows  one 
to  argue  kinematically  that  wake  waves  of  a  given  wavelength  are 
associated  with  the  SAR  signature  of  the  wake  because  they  have  the 
appropriate  group  propagation  velocity  and  were  presumably  generated 
near  the  ship.  For  the  interior  of  the  wake  near  the  ship  track  a 
small  change  in  wake  geometry  can  change  the  physics  one  employs  to 
explain  the  wake  image.  For  example,  a  change  from  2°  to  3°  in  the 
angle  0  ,  between  wake  and  ship  track,  changes  the  associated 
wavelength  by  more  than  a  factor  tw  -  from  20  era  to  45  cm.  Hence 
rather  complex  corrections  to  SAR  images  may  be  necessary — 
corrections  which  are  not  ordinarily  made.  During  June  1984,  Dr. 
David  Lysenga  and  Mr.  James  Lyden  of  the  Environmental  Research 
Institute  of  Michigan  (ERIM)  described  some  of  these  corrections  and 


their  applications  to  images  collected  during  the  Georgia  Straits 
Experiment.  Below  we  review  these  corrections  and  comment  upon  them. 

2.0  Scanning  Distortion 

This  correction  refers  to  the  finite  time  needed  to  collect 

the  SAR  image  of  a  ship  wake  and  the  fact  that  the  wake  waves  move 

during  this  time.  If  the  aircraft  flies  parallel  to  the  ship  track, 

the  correction  can  be  substantial.  It  requires  80s  for  a  200  kt 

aircraft  to  scan  a  10  km  wake  image.  During  this  time  radar  resonant 

(~30  cm)  ocean  waves  travel  some  56  m  or  about  6  SAR  resolution 

cells.  This  would  change  the  observed  wake  angle  0  by  about  0.2°. 

The  effect  is  illustrated  in  Figure  D-l  showing  that  the  scanning 

distortion  increases  the  observed  0  fov  SAR  movement  opposite  to  the 

ship  direction  and  vice-versa.  The  change  in  0  caused  by  this 

effect  is  0  ,  -  0..  =>  A0  «  -  c  /V  where  c„  is  the  group 

obs  true  g  8  43  r 

velocity  of  the  wake  waves  associated  with  the  SAR  signature  and  V  is 

the  aircraft  velocity  taken  positive  when  in  the  same  direction  as 

the  ship  velocity.  Clearly  scanning  distortion  is  much  less 

important  for  satellite  SAR  observations  since  V  ~  7,000  m/s  and 

0 

-3 

A0  for  ~30  cm  waves  becomes  only  ~3  x  10  • 
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3.0  Doppler  Shift 


Since  SAR  uses  the  Doppler  shift  of  the  scatterers  in  a 
resolution  cell  to  locate  that  resolution  cell  along  the  azimuth 
direction,  any  unexpected  radial  velocity  (unrelated  to  the  known 
movement  of  the  SAR  platform  with  respect  to  the  ocean  surface)  can 
cause  the  resolution  cell  to  be  shifted  along  the  azimuth  direction 
by  an  amount  Ay  =  (R/V)  where  R  is  the  range  to  the  resolution 
cell,  V  is  the  SAR  velocity  and  Vr  is  the  unexpected  radial 
velocity.  Bragg  resonant  waves  of  ~30°  have  a  phase  velocity  of 
M).?  m/s  and  thus  cause  a  typical  displacement  of 
Ay  ~(100)  0.35  ~35  m  at  an  incidence  angle  of  ~30°  .  If  a  ship 
wake  were  all  observed  at  a  single  angle  of  incidence,  the  wake  arms 
v'ould  simply  be  shifted  as  a  whole,  the  angle  between  the  arms  being 
unaffected.  Here  we  neglect  the  velocity  due  to  waves  in  the 
transverse  wake  since  their  velocity  is  primarily  perpendicular  to 
the  radar  ray  path  for  ships  traveling  nearly  along  the  SAR  flight 
direction.  However,  different  portions  of  the  wake  are  observed  at 
different  incidence  angles  and  therefore  interact  with  different 
wavelength  waves  which  in  turn  have  different  phase  velocities  and 
thus  different  Doppler  displacements.  Thus  a  small  distortion  of  the 
wake  angle  is  possible.  Since  Bragg  resonant  waves  at  larger 
incidence  angles  are  shorter,  they  travel  more  slowly;  yet  they  are 
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viewed  more  nearly  along  their  velocity  vector  and  have  a  larger 
component  along  the  radar  ray  path.  Since  these  two  effects  tend  to 
cancel  and  wakes  are  most  evident  at  small  incidence  angles  and 
nearly  parallel  to  the  SAR  velocity,  this  effect  is  usually  a  small 
distortion  compared  with  the  scanning  distortion  above.  The 
distortion  for  satellite  SAR  is  less  than  for  aircraft  SAR  since  the 
satellite  SAR  views  the  whole  ship  wake  at  very  nearly  the  same  angle 
of  incidence. 

4.0  Slant  to  Ground  Range  Conversion 

The  most  primitive  form  of  SAR  image  is  the  slant  range 
image  in  which  the  image  is  displayed  with  true  length  along  the  SAR 
flight  (y)  direction  but  with  the  other  dimension  being  range  (R) 
along  the  radar  ray  path  not  along  the  Earth's  surface.  The 
conversion  between  slant  and  ground  range  Is  simply  ground  range 
x  «  R  sin  0"  where  R  is  slant  range  and  0"  angle  of  incidence. 
Linear  features  which  are  straight  lines  In  the  ground  range,  (x,y) 
plane,  are  transformed  into  curved  lines  in  the  slant  range  (R,y) 
plane.  This  effect  is  illustrated  in  Figure  D-2.  For  aircraft  the 
distortion  can  be  striking  (especially  at  small  0'  because 
0"  usually  varies  over  a  large  range,  e.g.,  0°  to  55°  for  the 
aircraft  used  in  the  Georgia  Strait  experiment.  For  spacecraft  the 
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Figure  D-2.  Geometrical  distortion  in  the  slant-range  imagery 
(adapted  from  Fu  and  Holt,  1982). 
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effect  is  much  less  pronounced  because  only  a  small  range  of  6  is 
used,  e.g.,  19°  to  25°  for  SEASAT.  Aircraft  data  are  often  presented 
in  slant  range  images  since  they  are  more  easily  obtained.  Satellite 
data  on  the  other  hand  are  nearly  always  provided  in  ground  range 
coordinates.  Obviously  ground  range  images  must  be  used  to  measure 
wake  geometry. 

The  original  measurement  measurements  of  ship  wake  angles 
(8)  during  the  Georgia  Straits  (1983)  experiment  were  made  from 
slant  range  images  and  corrected  approximately  for  the  slant  to 
ground  range  conversion.  When  measurements  of  0  on  ground  range 
images  were  made  later,  it  was  found  that  the  original  approximate 
correction  used  with  slant  range  images — could  be  in  error  by  about 
0.7°.  This  effect  coupled  with  the  -^).2°  scanning  error  comprise 
the  ~1°  total  correction  overall. 

5.0  Variation  of  Radar  Cross  Section  with  Angle  of  Incidence 


The  radar  cross  section  o°  of  an  element  of  ocean  surface 
is  a  function,  not  only  of  the  local  surface  roughness,  but  also  of 
the  radar  parameters  and  observational  geometry.  In  particular 
varies  over  some  six  orders  of  magnitude  as  angle  of  Incidence 
varies  from  0°  to  near  90°.  So  if  the  ship  wake  being  observed  runs 
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over  a  range  of  0"  ,  some  distortion  could  be  introduced  because  of 
the  0(0')  variation.  That  is,  a  linear  wake  feature  would  be 
shifted  toward  lower  0"  since  da/d0  is  negative.  The  rate  of 
change  of  o  with  change  in  range  is 

do/dx  =  (da/d0')(d0"/dx)  =  (da/d0')  R_i  sec  Q'  .  For  a  typical 

aircraft  R  "10  km  and  da/d0"  “60  dB/rad  at  maximum 

( 0"  “20°),  do/dx  “  60  x  10~4  x  1.06  "  6  x  lO-3  dB/m  .  This  is  quite 

small,  <0.1  dB  for  a  one  resolution  cell  displacement.  So  we 

conclude  that  this  correction  is  negligible  for  aircraft  and  doubly 

so  for  spacecraft  (where  R  is  larger)  provided  0"  is  not  near  90°. 

0 

For  the  SAR's  of  interest  here  0  is  limited  to  <60  . 

6.0  Applications  to  Georgia  Strait  Experiment 

From  the  discussion  above  it  appears  that  the  most  important 
correction  is  to  put  the  image  in  ground  range  coordinates.  Next  in 
importance  is  the  scanning  correction  which  may  account  for  several 
tenths  of  a  degree.  The  other  corrections  appear  to  be  typically 
smaller.  David  Lyzenga  &  James  Lyden  of  ERIM  applied  the 
aforementioned  corrections  to  SAR  images  collected  during  the  Georgia 
Straits  Experiment.  These  corrections  caused  the  a r-s  of  the  wake  in 
pass  9-2  to  meet  some  400m  ahead  of  the  ship  image  as  one  would 
expect  from  Figure  D-l.  Corrections  to  narrow  V  wakes  observed 
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APPENDIX  E 


SPECULAR  POINT  SCATTERING  FROM  COHERENT  WAVE  PACKETS 

1.0  Introduction 

This  and  other  reports  (Zachariasen,  et  al.,  1984;  Hammond, 
et  al.,  1984,  and  Shuchman  et  al.,  1983)  describe  SAR  images  of  ship 
wakes.  These  reports  and  SAR  data  gathered  during  the  Gulf  of  Alaska 
SAR  -  ship  wake  experiment  have  revealed  a  number  of  characteristics 
common  to  narrow  angle  wakes  (6  ~  2.5°  in  Figure  E-l).  Some  of 
these  features  are  as  follows: 

1.  Location:  Wake  arms  are  located  at  Q  ~  2.5°  which 
corresponds  and  kinetically  to  ~  30  cm  waves 
propagating  outward  from  near  the  ship  track. 

2.  Length:  Wakes  are  very  long  extending  in  some  cases 
more  than  20  km  behind  the  ship. 

3*  Fineness :  The  arms  making  up  the  bright  wake  are  very 
fine,  occupying  a  few  SAR  resolution  cells  or  less, 

4.  Brightness :  The  bright  arms  making  up  the  wake  are 
indeed  very  bright,  being  like  6  to  24  dB  above  the 
surrounding  background  at  ~  1  km  behind  the  ship. 
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Figure  E-l.  Characteristics  common  to  narrow  angle  wakes 


Sensitivity  to  angle  of  Incidence:  Narrow  wakes  are 


very  prominent  at  small  incidence  angles  6  ~  20  >  but 

have  not  been  observed  at  high  incidence  angles 
0  ~  50°  •  Observation  of  the  same  ship  on  two 
successive  SAR  passes  a  few  minutes  apart 
at  0  ~  25°  and  0  ~  50°  is  convincing  evidence  for 

this  characteristic. 


Disappearance  in  High  Winds:  Narrow  angle,  SAR  ship 
wakes  are  most  prominent  at  very  low  winds  and  are 
either  very  much  reduced  in  length  or  disappear  in  high 
winds . 


In  this  section  we  discuss  an  hypothesis  for  the  scattering 
mechanism  which  addresses  how  the  above  characteristics  of  ship  wakes 
occur  in  SAR  images.  The  essence  of  the  hypothesis  is  that  the  radar 
energy  is  backscatted  from  specular  points  located  in  small  regions 
of  coherent  ~  30  cm  (or  possibly  ~  60  cm)  surface  waves.  These  small 
regions  are  generated  as  part  of  the  ship's  surface  (Kelvin)  wake. 
Since  the  waves  involved  are  short  and  of  height  <  few  centimeters, 
they  are  not  particularly  noticeable  in  a  photograph  of  the  wake. 
However,  because  microwave  radar  is  sensitive  to  length  scales  of  cm 


to  10's  of  cm  and  particularly  so  to  coherent'  patches  of  such  waves, 
we  suggest  that  these  waves  form  the  narrow  vee  SAR  signature  of  a 
ship  wake. 

2.0  SPECULAR  POINT  SCATTERING  FROM  COHERENT  WAVE  PACKETS 

2.1  Coherence  of  Wake  Waves  Within  A  Resolution  Cell 

One  of  the  most  striking  visual  features  of  a  ship  wake  in 
calm  water  is  the  coherence  of  waves  in  the  transverse  and  divergent 
portions  of  the  wake.  The  divergent  portion  of  the  wake  is  of 
interest  here  since  the  SAR  images  of  narrow  wakes  are  for  wake  arms 
which  are  nearly  parallel  to  the  SAR  velocity  vector,  i.e.,  0 
and  6  -  a  few  degrees  in  Figure  E-l.  The  cover  photo  of  USS  Quapav 
(Hammond,  et  al.,  1984)  demonstrates  this.  Using  the  classical 
Kelvin  wake  analysis  summarized  by  Dashen  et  al.  (1982).  We  can 
calculate  the  variation  in  wave  number  between  given 
angles  0  »  0^  and  0  ■  0^  in  Figure  E-l.  In  Figure  E-2  we  show  the 
variation  in  wavelength  with  0  over  several  ranges  of  0. 


Wake  Vertex  Angles 
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Figure  E-2. 


Wake  Vertex  Angles 


Variation  of  classical  Kelvin  wavelength  with  angle  0  (between  ship  track  nnd 
line  at  angle  to  ship  track  -  0  -  ±19.4*  shows  cusp  lines  of  classical 
Kelvin  wake).  In  (a)  0  ranges  from  0°  to  5a.  In  (b)  0  ranges  from  2°  to  3° 
and  we  show  the  range  of  angles  and  hence  the  range  of  wake  wavelengths 
within  a  10m  SAR  resolution  cell  at  1,  2,  5  and  10km  distances  behind  the 


range  of  wavelengths  (A)  around  the  31  cm  wavelength  at  the 
8  «  2.5°  point  is  -7  cm  (0.23A),  -3.5  cm  (0.11A), 

-1,5  cm  (0.05A),  and  -0.75  cm  (0.02A)  at  distances  of  1,  2,  5,  and 
10  km  behind  the  ship  (x  in  Figure  E-l).  The  wavelength  within  a  SAR 
resolution  cell  becomes  more  pure  at  greater  distances  behind  the 
ship  because  the  resolution  cell  subtends  smaller  and  smaller 
variations  of  9.  Remember,  of  course,  that  these  calculations  refer 
to  an  ideal  classical  Kelvin  wake  in  which  the  ship  is  a  point 
source. 


It  is  common  to  characterize  the  sea  surface  height 
fluctuations  by  a  two-dimensional  Fourier  power  spectrum  (Phillips, 
1977).  Such  a  characterization  assumes  that  the  ocean  surface  waves 
are  randomly  phased  as  indeed  they  are  for  a  wind  roughened  sea.  In 
the  case  of  ship  wakes  the  waves  are  clearly  not  randomly  phased  in 
general,  but  rather  coherent.  Since  the  customary  Fourier  power 
spectrum  throws  away  phase,  it  can  not  adequately  characterize  the 
waves  in  a  ship  wake.  Other  statistics  which  can  characterize  a 
partially  coherent  phenomena  are  necessary. 

Perusal  of  the  cover  photo  and  Figure  19  of  Hammoud  et  al. 
(1934)  quickly  shows  that  for  a  real  wake  at  0  «■  2.5°  there  are 
substantial  deviations  from  the  classical  Kelvin  wake  picture. 


First,  the  wake  waves  have  a  significant  distribution  in  direction 
appearing  to  radiate  from  a  sequence  of  compact  sources  spaced 
approximately  at  the  wake's  transverse  wavelength.  Thus  the  radius 
of  curvature  (in  the  horizontal  plane)  appears  not  to  be  that  of  the 
divergent  portion  of  the  classical  Kelvin  wake,  but  rather  ~  the 
distance  from  the  ship  track  to  the  point  of  interest.  Second,  the 
dominant  wavelengths  at  a  given  location  deviate  from  the  classical 
Kelvin  wake  predictions  of  Figure  E-2.  Finally,  the  waves  appear 
only  partially  coherent.  The  characteristics  of  wake  waves  at 
small  0  suggested  by  Figure  19  are  consistent  with  the  intermittent 
source  picture  discussed  in  section  3  of  this  report. 

Although  this  picture,  suggested  by  observations,  is 
somewhat  different  from  the  idealized  classical  Kelvin  wake;  features 
important  for  explaining  the  SAR  signature  of  ship  wakes  are 
present.  In  fact  some  observational  features  of  ship  wake  waves  are 
necessary  to  explain  the  SAR  observations.  For  example,  the 
directional  spread  of  ~  30  cm  is  necessary  to  allow  both  arms  of  a 


narrow  vee  wake  to  be  observed 


2 . 2  Radar  Cross  Section  From  Coherent  Specular  Scatter ers 


1 .  Definition  of  radar  cross  section 

For  people  unfamiliar  with  radar,  the  notion  of  radar  cross 
has  some  unexpected  features  when  specular  scattering  is  present. 

The  definition  of  radar  cross  section  a  is  given  by 

-  Power  reflected  toward  observer/unit  solid  angle 
incident  power  density/4  ir 


where  R  is  the  distance  between  the  radar  and  the  scatterer,  Er  is 
the  reflected  field  strength  at  the  observer  and  ^  is  the  incident 
field  strength  at  the  scatterer.  The  radar  cross  section  is  the  area 
necessary  to  intercept  the  amount  ot  incident  power  which,  when 
scattered  isotropically,  produces  the  echo  power  actually  observed. 
Thus  a  is  not_  a  real  cross  sectional  area,  but  rather  a  convenient 
hypothetical  area.  When  o  is  measured  according  to  equation  (E-l)  it 
can  be  used  in  equation  (E-2)  to  predict  the  backscatter  power  from  a 
target,  namely 


p  .w 

r  (4*)2  R* 


(E-2) 


E-8 


where  Pr  and  Pt  are  the  received  and  transmitted  power,  Gt  is  the 
transmitter  antenna  gain,  A£  is  the  effective  area  of  the  receiver 
antenna,  R  is  the  range  and  a  the  radar  cross  section  (in  units  of 
area) .  When  land  or  water  surfaces  are  being  observed  the 
quantity  o°  denotes  the  radar  cross  section  a  of  a  unit  area, 
i.e.,  o°  =  o/A  where  a  is  the  cross  section  from  area  A. 

2.  Radar  cross  section  of  specular  reflectors: 

The  counter  intuitive  property  of  o°  is  that  it  can  be  and 

often  is  greater  than  unity.  For  example  a  lossless  reflecting 

surface  of  area  A  oriented  for  specular  reflection  has  a  cross 

2  2 

section  of  o  ~  (  4ttA  /X  ).  Note  that  a  is  proportional  to  the 
physical  area  of  the  reflector  squared  because  the  scattering  is 
coherent  over  the  surface. 

3 *  Kodis  model  for  radar  cross  section  of  a  sinusoidal 
corrugated  surface 

Kodis  (1966)  discusses  the  radar  cross  section  from  surfaces 
containing  specular  points.  The  theory  is  applicable  to  surfaces 
where  the  "physical  optics"  or  "Kirchoff"  method  can  be  used,  i.e., 
where  pk  cos3  G  »  1  (E-3) 

(Brehkovskikh,  1952).  Such  a  surface  is  called  "gently  undulating." 


E-9 


Consider  a  sinusoidally  corrugated  surface  of  width  W  along 
the  crest  of  a  corrugation  and  length  L  perpendicular  to  the 
corrugation  crests.  The  radar  illuminates  this  corrugated  surface  in 
a  direction  perpendicular  to  the  corrugation  crests  so  that  specular 
backscatter  may  arise.  If  the  angle  of  incidence  0  is  sufficiently 
small  (near  normal  incidence)  the  the  slope  of  the  sinusoidal  surface 
provides  specular  points.  The  radar  cross  section  is  then  governed 
by  the  phasing  of  the  specular  points,  i.e.  for  backscatter  the  radar 
cross  section  per  unit  width  is  (after  Kodis) 

n-1 

a  *  (irp  /X)  [N  +  2  1  (N  -  m)  cos  m  <J>]  (E-4) 

8  n-1 

where  p  is  the  radius  of  curvature  of  the  sinusoidal  surface  at  the 
s 

specular  point,  X  is  the  radar  wavelength,  N  the  number  of  specular 
point  illuminated  and  <(>  =  2  kAn^  where  k  is  the  radar 
wavenumber,  A  the  sinusoidal  surface  wavelength  and  nx  is  the 
component  of  the  surface  normal  along  the  illumination  or  x- 
direction.  As  with  a  grating,  backscatter  maximum  occur  when  <{>  is  an 
integral  multiple  of  2v  leading  to  a  Bragg  condition 


When  the  Bragg  condition  is  fulfilled,  the  radar  cross  section  for 
backscatter  from  an  area  of  length  L  and  width  W  becomes  a  maximum, 
namely  (after  Kodis) 

0  “  <2*PaA)  N2  (LW)  (E-6) 

max  s 

where  N  is  the  number  of  specular  points  illuminated . 

For  a  randomly  rough  surface  fulfilling  (3)  Kodis  obtains  a 
backscatter  cross  section  of 

0  ~  t r  R~R^  <N>  (E-7  ) 

where  Ra  and  R^  are  the  principle  radii  of  curvature  of  a  specular 
point  and  <N>  is  the  mean  number  of  specular  points  in  the 
illuminated  region,  a  becomes  0°  if  we  take  <K>  to  be  the  mean 
number  of  specular  points  per  unit  area.  The  important  point  to  note 
here  is  that  o  goes  as  for  the  coherent  corrugated  surface,  but  as 
N  for  the  randomly  rough  surface . 


2.3  Radar  Cross  Section  for  a  SAR  Resolution  Cell  Containing 
Coherent  Specular  Scatterers 

We  now  apply  the  notions  above  to  a  SAR  resolution  cell 
containing  at  least  some  portion  which  is  filled  with  coherent  waves 
as  discussed  above.  We  consider  sinusoidal  waves  in  some  detail  and 
discuss  other  wave  shapes  briefly. 

1.  Sinusoidal  Waves. 

Kodis  (1966)  Consider  a  sinusoidal  corrugated  surface 

of  the  form 

C(x)  =  h  sin  (2itx/&)  (E-8 

where  — (L/2)  <  x  <  (L/2).  The  maximum  slope  of  the  surface 
is  (2ir  h/£)  so  that  specular  points  will  appear  only  for  angles  of 
incidence  Q  smaller  than 

tan  0  ■  2nh/ %  .  (E-9 

The  radius  of  curvature  at  the  specular  points  is 


P„  a 


A^/(4n  h  cos^  0  sin  <j;) 


(E-10 


where  cos  4)  =  (Jt/2n  h)  tan  8  .  In  general  there  are  twc  specular 

points  such  that  (2ir/£)  x  =  \|>  and  (2nM)  x  „  =  2ir  -  <J).  The  cross 

si  s2 

section  of  equation  (E— 6)  assumes  that  only  one  of  the  two  is 
illuminated,  i.e.,  one  is  shadowed  or  the  two  coincide  at  the  poiut 
of  maximum  slope. 

As  a  specific  case,  we  take  8  =  23°  as  with  SEASAT 

and  X  =  23  cm  for  an  L-band  SAR  yielding  a  Bragg  resonant 

Z  of  ~  30  cm.  In  Table  E-l  we  note  the  minimum  wave  amplitude  (h  ) 

s 

required  to  produce  a  specular  point  at  the  point  of  maximum  slope  on 
the  wave.  Also  shown  is  the  maximum  slope  for  the  required  h8 
and  Z  ~  30  cm.  Since  waves  tend  to  break  when  they  have  maximum 
slopes  above  ~  0.5,  it  is  clear  that  for  sinusoidal  waves,  angles  of 

A 

incidence  8  above  about  30°  will  not  produce  specular  points  and 

•fa 

thus  no  radar  backscatter  via  this  mechanism.  We  also  show  p  for 

s 

conditions  such  that  4»  =  180°  -  150°,  i.e.,  near  but  not  at  the  point 

of  maximum  slope  (4)  =  0°).  If  the  specular  points  are  precisely  at 

the  point  of  maximum  slope,  then  sin  \Jj  =  0  and  p  becomes  infinitely 

s 

large. 

In  a  10  x  10  ra  resolution  cell  with  p  ~  0.9  m  ,  and 

s 

X  ~  0.23  m  we  have  from  equation  (E-6'<  and  o°  =>  o/Area 


o 


a 
max 


(E-ll) 


where  N  is  the  number  of  coherent  specular  points  matching  the  Bragg 
condition  (equation  (E-5)).  Referring  to  Figure  E-3  we  see  that  even 
a  few  coherent  corregations  will  produce  a  normalized  radar  cross 
section  o°  which  is  greater  than  the  observed  radar  cross  section  of 
the  ocean  by  at  least  10  dB. 


Now  equation  (E-6)  refers  to  a  straight  corrugation 
stretching  clear  across  the  resolution  cell.  As  we  noted  above,  the 
observed  ~  30  cm  waves,  have  significant  curvature  in  the  horizontal 
plane  and  clearly  do  not  have  wave  crests  10  m  long.  Further,  one 
might  ask,  is  the  illumination  coherent  over  the  resolution  cell? 
Zachariasen  et  al.  (1983)  note  that  for  all  common  SAR  systems  the 
Fresnel  zone  size  is  larger  than  the  resolution  cell  size.  If  we 
accept  the  idea  that  the  ~  30  cm  waves  are  generated  by  an 
intermittent  source  near  the  ship  track,  then  the  radius  of  curvature 
in  the  horizontal  plane  must  be  about  the  distance  from 
the  ship  track  to  the  wake's  a  SAR  signature  at  9  ~  2.5° 
i.e.  R^  "  0.04x  or  R^  ~  44  m  at  one  kilometer  behind  the  ship. 

The  coherent  length  of  wave  crest  is  l  ~  or  if  X  ~  0.23  m 

and  Rj.  ~  44  m,  ~  2.2  m.  Thus  the  o°  from  equation  (E-ll)  should 


E— 15 


be  lowered  at  least  by  the  factor  (£  /W)  where  W  is  the  resolution 

c 

cell  width.  Thus  equation  (2.11)  becomes: 


o°  ~  5.5  N2 
max 


(E-12) 


which  is  still  more  than  10  dB  above  the  level  of  Figure  E-3. 


In  summary  the  point  of  this  computation  is  to  show  that 
even  a  small  patch  of  several  coherent  waves  with  specular  points  can 
provide  a  radar  cross  section  significantly  above  that  of  the  average 
ocean , 


2.  Cynoidal,  Gusped  and  Breaking  Waves. 

Realistic  surface  waves  can  have  cynoidal  or  cusped 
shapes.  These  shapes  have  regions  with  radii  of  curvature  smaller 
than  the  radar  wavelength  and  thus  can  be  modeled  approximately  as 
scattering  wedges.  Some  work  on  wedge  scattering  in  connection  with 
the  ocean  has  been  done  by  Lyzenga  et  al  (1983)  and  a  laboratory  wave 
tank  by  Lake  and  Kwoh  at  TRW.  Waves  of  ~  30  cm  lengths  with  these 
shapes  could  also  produce  specular  scattering  from  coherent  wave 
packets  either  via  specular  points  or  by  wedge  shaped  crests.  The 
radar  cross  section  from  these  shapes  would  clearly  differ  from 
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equations  (E-ll  or  E-12),  but  if  the  waves  were  in  coherent  sets  the 


total  cross  section  would  go  as  rather  than  N.  Hence  sets  of 
coherent  wedge  scathers  could  produce  o°  significantly  above  that  of 
the  average  ocean. 

3.  Importance  of  Polarization 

Nearly  all  SAR's  observing  the  ocean,  such  as  SEASAT 
and  the  NASA  L-band  aircraft  SAR,  use  HH  polarization.  Thus 
horizontally  polarized  (E  field  perpendicular  to  the  plane  of 
incidence)  waves  are  transmitted  and  received.  Since  the  coherent 
waves  discussed  in  this  section  have  much  longer  coherence  lengths  in 
the  horizontal  direction,  one  would  expect  a  stronger  wake  signature 
in  HH  as  opposed  to  VV  polarization.  This  effect  would  be 
particularly  striking  if  wedge  like  structures  were  present.  So  SAR 
images  of  ship  wakes  simultaneously  in  both  HH  and  W  (and  possibly 
HV  or  VH  polarization)  would  carry  important  information  regarding 
the  wave  structures  responsible  for  the  observed  scattering. 

2.4  Validity  of  Coherent  Specular  Model 

The  key  question  concerning  the  validity  of  the  hypothesis 
suggested  here  is  what  are  the  wake  waves  responsible  for  the  SAR 
signature  really  like?  Do  they  satisfy  the  curvature  condition 
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equation  (E-3)?  Do  they  contain  specular  points?  Are  wedge-like 
structures  involved?  How  coherent  are  the  real  waves  in  nature? 


This  coherent  specular  hypothesis  rests  on  the  conjecture  that  real 
wake  waves  at  ~  30  cm  length  are  enough  like  an  ideal  coherent- 
sinusoidally-corrugated  surface  that  a  large  radar  cross  section  is 
produced . 

3.0  APPLICATION  TO  SAR  WAKE  FEATURES 


Below  we  discuss  one  by  one  the  observed  features  of  ship 
wakes  observed  by  L-band  SAR.  In  each  case  we  show  how  this 
coherent-specular-scatter  dypothesis  can  explain  or  help  explain  the 
observed  signature.  Since  the  scattering  surface  is  not  well 
characterized,  these  are  arguements  are  at  best  semi-quantitative. 

3.1  Location  of  Wake 


SAR  images  at  L-band  find  bright  wake  signatures 
at  0  ~  2.5°  (Figure  E-l).  This  value  of  0  corresponds  to  the  time 
required  for  ~  30  cm  waves  to  propagate  from  the  ship  track 
to  0  ~  2.5°.  Thus  one  can  argue  kinematically  that  short  gravity 
waves  which  fulfill  the  Bragg  resonance  condition  of  equation  (E-5) 
are  present  at  the  proper  location.  Hammond  et  al.  (1984)  find  from 
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photographic  data  that  the  wave  energy  present  at  6  ~  2.5°  is 
concentrated  more  at  wavelengths  of  from  50  to  90  cm  than  from  10  to 
50  cm.  This  results  is,  not  comprehensive,  however,  being  composed 
of  one  observation.  Nevertheless  60  cm  waves  also  fulfill  the  Bragg 
condition  of  equation  (E-5).  They  would  produce  a  lower  (but  still 
large)  cross  section  because  N  would  be  smaller  in  equations  (E-ll) 
and  CE— 12) . 


3.2  Length  of  Wake 

SAR  images  of  ship  wakes  stretch  to  20  km  or  more  along  the 
x-direction  in  Figure  E-l.  Given  that  30  cm  waves  decay  with  an  e- 
folding  time  considerably  less  than  the  ~  30  min  required  to  generate 
such  a  long  wake.  One  suspects  that  there  must  be  some  mechanism 
which  enhances  the  radar  cross  section  to  at  least  partially 
counteract  the  wave  decay.  The  coherent  specular  hypothesis  provides 
such  a  mechanism.  In  Figure  E-2  we  showed  how  the  spread  in 
wavelength  of  classical  Kelvin  wake  waves  within  a  SAR  resolution 
cell  becomes  smaller  as  the  resolution  cell  in  question  was  located 
further  behind  the  ship.  Thus  as  x  increases  and  **  30  cm  waves 
decay,  the  purity  of  the  wavelength  within  a  10  x  10  m  resolution 
cell  increases.  Hence  the  fraction  of  the  resolution  cell  covered  by 
a  resonant  corrugated  surface  would  increase,  N  would  Increase 


and  0  would  increase  according  to  equation  (E-6).  However, 

as  ~  30  cm  waves  decay,  they  become  less  able  to  fulfill  equation 

(E-9)  and  thus  produce  specular  points  for  a  given  angle  of 
* 

incidence  0  . 

3.3  Fine  (Narrow)  Wake  Features 


The  fact  that  SAR  ship  wake  features  occur  only  one  to  a  few 
SAR  resolution  cells  suggests  that  a  sharp  resonance  phenomena  is 
present.  The  coherent  specular  hypothesis  fits  this  observational 
fact  since  the  Bragg  condition  (equation  (E-5))  must  be  fulfilled  in 
order  to  obtain  the  maximum  cross  section  of  equation  (E-6). 
Scattering  via  the  slightly  rough  Bragg  resonant  mechanism 
(Valenzuela,  1979)  also  fits  well  with  the  fineness  of  the  linear 
wake  features.  The  coherent  specular  mechanism  differs  from  the 
slightly  rough  mechanism  in  that  the  former  depends  on  specular 
scattering  from  coherently  phased  waves  while  the  later  does  not 
depend  on  specualr  scattering.  An  analysis  of  the  slightly  rough 
Bragg  mechanism  needs  to  be  done  to  compare  it  with  the  hypothesis 
pursued  here.  However,  such  a  comparison  is  beyond  the  scope  of  this 


report 


3.4  Strong  Wake  Signature 


The  observation  that  SAR  wake  signature  in  calm  water  are 
very  strong  (~  8  to  18  dB  above  the  background)  suggests  that  a 
strong  scattering  mechanism  is  required.  In  section  2  above  we 
showed  how  even  a  small  N  in  equation  (E-12)  could  provide  a  radar 
cross  section  some  10  to  20  dB  above  the  normal  ocean  cross  section 
of  Figure  E-3.  In  F  below  we  discuss  why  such  a  large  enhancement 
in  o  is  not  seen  when  the  wind  rises  and  other  waves  beside  wake 
generated  waves  are  present. 

3.5  Disappearance  of  Narrow  Vee  SAR  WAkes  at  High  Incidence  Angles 


During  the  Gulf  of  Alaska  SAR  ship  wake  experiment  the  same 
ship  was  observed  at  o  ~  50°  the  same  ship  in  the  same  background 
produced  no  wake  in  the  L-band  SAR  image.  The  coherent  specular 
hypothesis  explains  this  qualitatively  in  Table  E-l.  At 
large  a  large  wave  slopes  are  required  to  produce  specular  points. 
For  a  ^  30°  the  Bragg  resonant  waves  break  before  they  can  produces 
the  large  required  slopes.  Further,  as  a '  increases  the  Bragg 
resonance  condition  requires  shorter  surface  waves  for  resonance. 

At  o  «  5C°  the  first  order  resonance  is  at  15  cm  rather  than  30 
cm.  The  shorter  waves  decay  faster  and  may  be  generated  at  less 
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strength  and  less  coherently  by  the  ship  -  more  reasons  why  narrow 
SAR  wakes  are  not  observed  at  large  <j  « 


3 .6  Disappearance  of  Narrow  Wakes  as  Wind  Speed  Increases 

If  coherent  specular  scattering  can  produce  such  a 
large  o  on  a  calm  surface,  why  can't  this  mechanism  make  narrow  wakes 
visible  in  even  high  sea  states.  We  argue  that  as  the  wind  increases 
not  only  does  the  background  a  (surrounding  the  location  of  the 
narrow  vee  wake)  rise  but  also  that  the  wind  generated  waves  distroy 
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